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20. . Abstract  (continued) 

A detailed  model  of  a sunspot  penumbra  is  studied  in  order  to  identify 
the  mode  of  running  penumbral  waves.  --Mt  is  found  that  penumbral  waves  may 
be  identified  with  magneto-atmospheric  waves  of  the^pl us^type  that  are 
vertically  trapped  at  photospheric  levels.  Although  most  of  the  wave  energy 
is  contained  in  the  penumbral  photosphere  and  subphotosphere,  the  maximum 
vertical  velocity  occurs  in  the  chromosphere  where  the  waves  are  evanescent 
(and  where  they  are  observed  in  Ho). 


An  exact  analytical  solution  for  magneto-atmospheric  wave  modes  is  found 
in  the  case  of  an  isothermal  atmosphere  permeated  by  a uniform  horizontal 
magnetic  field,  without  making  the  usual  short-wavelength  approximation. 

This  solution  is  applied  to  an  idealized  model  of  the  low-corona-chromosphere 
transition  region  as  a model  for  flare-induced  coronal  waves.  Disturbances 
propagate  horizontally  in  the  waveguide  formed  by  the  rapid  density  increase 
into  the  chromosphere  below  and  the  rapid  increase  in  Alfv6n  speed  into  the 
corona  above. 


'-‘‘The  exact  solution  mentioned  above  is  also  used  in  conjunction  with  a 
simple  two-layer  model  of  a sunspot  penumbra  to  further  study  the  mode  of 
running  penumbral  waves.  The  lowest -'plus*'*  eigenmode  of  the  model  is  in 
good  agreement  with  observations  of  penumbral  waves.  ^ 


The  theory  of  penumbral  waves  developed  here  predicts  that  these  waves 
should  be  observable  in  the  photosphere  as  well  as  the  chromosphere.  This 
prediction  prompted  a search  for  penumbral  waves  in  the  photosphere,  carried 
out  with  the  tower  telescope  and  diode  array  at  Sacramento  Peak  Observatory. 
Simultaneous  observations  have  been  made  of  velocities  in  the  chromosphere 
(in  Ha)  and  in  the  photosphere  (in  the  nonmagnetic  Fe  I line  X5576)  of 
three  sunspots.  The  results  reveal  waves  propagating  horizontally  outward 
across  the  penumbra  in  the  photosphere  with  about  the  same  period  as  the 
running  penumbral  waves  in  Ha  (250-290  s).  The  photospheric  waves  are  more 
intermittent  and  have  higher  horizontal  phase  velocity  (by  a factor  of  2 or 
more)  than  the  chromospheric  penumbral  waves.  The  connection  between  the 
photospheric  and  chromospheric  penumbral  waves  is  unclear  at  present,  and  in 
any  case  is  more  complicated  than  the  resonant  mode  model  presented  here. 
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THE  NATURE  OK  RUNNING  PENUMBRAL  WAVES 


ALAN  H.  NYh*  and  JOHN  H.  THOMAS" 

Mus-PUmck-lnstitut  / in  Pltysik  mul  Astrophysik.  Munich,  Germany 

(Received  22  March;  revised  25  June,  1974) 

Abstract.  A model  of  a sunspot  penumbra,  including  (he  etTecis  of  magnetic  field,  compressibility,  and 
buoyancy,  is  studied  in  order  to  identify  the  mode  of  running  penumbral  waves.  It  is  found  that  the 
penumbral  waves  may  be  identified  with  gravity-modified  magneto-acoustic  waves  of  the  plus’  type 
that  arc  vertically  trapped  at  photospheric  levels.  Although  most  of  the  wave  energy  is  contained  in 
the  penumbral  photosphere  and  subphotosphere,  (he  maximum  vertical  velocity  occurs  m (he  chro- 
mosphere where  (i)  the  waves  arc  evanescent  and  (ii)  the  vertical  velocity  is  in  fact  observed  (in  Hu. 


1.  Introduction 

Recent  observations  have  disclosed  an  interesting  pattern  of  velocity  fields  in  sun- 
spots. The  most  recent  discovery  is  that  of  waves  propagating  radially  outward  in 
sunspot  penumbrae  (Zirin  and  Stein,  1972;  Giovanelli,  1972).  Zirin  and  Stein  refer  to 
these  wav  es  as  running  penumbral  waves.  With  the  further  observations  of  Giovanelli 
(1974),  we  now-  have  a fairly  clear  picture  of  the  properties  of  these  waves.  The  purpose 
of  this  theoretical  paper  is  to  study  possible  wave  modes  in  a model  of  a sunspot 
penumbra  in  order  to  identify  the  mode  of  the  running  penumbral  waves.  We  shall 
argue  that  the  running  penumbral  waves  are  gravity-modified  magneto-acoustic 
waves  (of  the  'plus’  type)  that  are  vertically  trapped  at  photospheric  levels. 

Giovanelli  (1974)  has  summarized  the  observations  of  running  penumbral  waves, 
and  he  presents  the  following  picture.  The  waves  are  observed  in  Hot  by  means  of 
their  line-of-sight  velocity.  They  occur  in  almost  every  sizable  spot  with  a regular 
stable  structure,  but  only  rarely  in  active  spots  with  complex  structure.  The  waves 
travel  outward  in  the  penumbra  at  a typical  speed  of  15  km  s"‘.  The  observed  waves 
have  periods  in  the  range  180  -240  s and  horizontal  wavelengths  in  the  range  2350- 
3800  km.  Observations  near  the  limb  have  failed  to  reveal  any  horizontal  motions 
associated  with  the  penumbral  waves,  so  the  wave  motion  is  predominantly  vertical 
in  Ha. 

Thus  far  no  detailed  theoretical  study  of  the  mode  of  the  running  penumbral  waves 
has  appeared,  although  Moore  (1973)  has  studied  the  related  problem  of  the  genera- 
tion of  penumbral  waves  in  the  umbra.  Zirin  and  Stein  ( 1972)  (tentatively  identified 
the  penumbral  waves  as  sound  waves,  whereas  Giovanelli  (1972,  1974)  identified  them 
as  Alfven  waves.  The  penumbral  waves,  with  their  predominant  vertical  motions,  no 
uoubt  involve  the  combined  effects  of  restoring  forces  due  to  compressibility,  magnetic 


* Also  National  Science  Foundation  Predoctoral  Trainee,  Dept,  of  Mechanical  and  Aerospace 
Sciences,  University  of  Rochester,  Rochester,  N.Y.,  U.S.A. 

**  On  eavc  of  absence  from  the  Dept,  of  Mechanical  and  Aerospace  Sciences  and  the  C.  E.  Kenneth 
Mees  Observatory,  University  of  Rochester,  Rochester,  N.Y.,  U.S.A. 
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field,  and  buoyancy,  and  a complete  theory  should  account  for  all  three  effects.  This 
is  done  in  the  present  paper. 

In  studying  penumhral  waves,  we  face  a difficulty,  in  that  there  seems  to  be  no 
complete,  generally  accepted  penuntbral  model  on  which  to  base  our  calculations. 
We  have  therefore  constructed  a penuntbral  model  for  use  in  studying  wave  modes. 
This  model  (presented  in  Section  3),  while  simple  enough  to  permit  analysis  of  wave 
modes,  nevertheless  reproduces  all  of  the  relevant  features  of  penuntbral  structure, 
and  is  in  reasonable  quantitative  agreement  with  observations.  We  have  assumed  the 
penumbral  magnetic  held  to  be  purely  horizontal,  but  varying  with  height.  True  pe- 
numbral  magnetic  fields  are  not  purely  horizontal,  although  they  may  be  very  nearly 
so  (Nishi  and  Makita,  1973).  There  is  some  disagreement  over  the  inclination  of  the 
magnetic  field  in  a penumbra  (see  Beckers  and  Schroter  (1969)  for  a summary  of 
observations).  The  assumption  of  a horizontal  field  here  is  mostly  a matter  of  con- 
venience; the  basic  mechanism  we  propose  for  the  vertical  trapping  of  penumbral 
waves  will  also  work  for  an  inclined  field.  We  have  also  taken  our  model  to  be  horizon- 
tally uniform  - that  is,  we  have  not  tried  to  represent  the  horizontal  filamentary 
structure  of  a penumbra  or  the  radial  geometry. 

In  Section  2 we  present  the  basic  equations  for  waves  in  our  model  penumbral 
atmosphere.  The  basic  atmosphere  is  completely  characterized  in  these  equations  by 
the  vertical  distribution  of  three  parameters:  the  sound  speed  c,  the  Alfven  velocity 
i\,  and  the  local  density  scale  height  H.  In  order  to  illustrate  the  properties  of  the 
various  wave  modes  that  can  occur,  we  study  the  dispersion  relation  that  holds  in  the 
case  of  constant  c,  oA,  and  H.  In  Section  3 we  present  the  basic  penumbral  model  in 
terms  of  the  distributions  of  c,  i\,  and  //  w ith  height.  In  Section  4 we  show  that  the 
penumbral  waves  may  be  identified  with  'plus'  modes  that  are  trapped  in  the  photo- 
spheric-subphotospheric  region  in  our  model.  We  discuss  these  modes  further  in 
Section  5. 

2.  Basic  Equations  and  Dispersion  Relations 


In  our  simplified  treatment  of  a sunspot  penumbra  we  shall  ignore  the  radial  spreading 
of  magnetic  field  lines,  and  consider  the  undisturbed  magnetic  field  to  be  purely 
horizontal  (in  the  jc-direction)  and  varying  with  height  z:  i.e.,  B„  = («„(.-),  0.0). 
We  assume  the  field  permeates  an  inviscid,  perfectly  conducting,  plane  siratilied 
atmosphere  with  constant  acceleration  of  gravity  g (=0.274  km  s'2)  in  the  negative 
e-direction.  The  undisturbed  pressure,  density,  and  temperature  are  denoted  by 
Po(z)<  <Jo(z )•  and  To(z)'  respectively.  The  atmosphere  is  in  hydrostatic  equilibrium, 
so  that 


d 


dz 


(1) 


We  then  consider  small  adiabatic  perturbations  of  this  equilibrium  atmosphere. 
We  consider  wave  vectors  only  in  the  .vz  plane,  and  assume  that  the  perturbation 
velocity  u =(u.v,w)  has  the  form  u = u exp/'  (kxx—(ot ),  with  u = u(z)=  (w(z),  f(z), 
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»■(;)).  Starting  with  the  linearized  equations  of  continuity,  energy,  and  momentum, 
we  can  eliminate  the  pressure  and  density  perturbations,  and  arrive  at  the  set  of  three 
linearized  momentum  equations  for  the  velocity  components  u,  P,  and  w1.  in  which  the 
basic  atmosphere  is  complete!)  described  In  the  sound  speed  i'(z),  the  Allven  velocity 
rA(r).  and  the  local  density  scale  height  II  (z),  delined  by 


The  linearized  momentum  equations  are  the  following: 


(w2  - c2k2x)  u 

, dw 

+ ikxc  — ikxgw  = 0, 

(3) 

dz 

(w2  - vik2)  v 

, , d2w 

{C  + & dz2 

= 0, 

Td  , , c2  + ofl 

+U c +”  J 

d/  + (o >2-vik2x)*  + 
dz 

(4) 

r dc2  c2 

1 , dii 

(5) 

+ ikx  — 

+ g w + ikxc2  = 0 . 

Ldz  H 

: dz 

We  can  eliminate  the  horizontal  velocity  components  ti  and  v from  the  system  of 
Equations  (3)— (5)  to  obtain  a single  equation  for  the  vertical  velocity  w,  in  the  form 


where 


and 


d2H’  dw 

2 + A(z)  — + B(z)w  = 0, 
dr  dr 


1 r co 

a<’\~h+  L(^ 

X [vl(w2-c2kl)  + cWyl 


dc2  , 2 2,2\dv 

l,s  , +(W2  -C  kx) 

c ki)  d z dz 


3- 


B( 


) (oj  - c2k2x) -g[g-  )kx- 
H / 


z)  = Qco2  - vik 

-9*  ? —ijj  ^.-1  0a  («2  - c2k2x)  + c2w2] 

(cr  - c kx)  d z J 


With  the  transformation 


(6) 


(7) 


(8) 


<P(z)  = w(z)exp  J zl(z)dzj 
0 

Equation  (6)  assumes  the  form 


d 2d> 

,2+J  (kx,  w;  z)<p  = 0, 
dz 


(9) 


(10) 
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where 

d/4 

f = B — \A‘  — | . (II) 

dz 

Here,  |<)>|2  is  roughly  proportional  to  the  kinetic  energy  of  wave  motion.  Equation 
(10)  is  the  propagation  equation  for  waves  in  an  atmosphere  with  vertically  varying 
c2,  vi,  and  II.  For  given  distributions  of  c2  (z),  u*  (z),  and  H (z),  we  can  use  the  expres- 
sion for /(z)  (Equation  (II))  to  distinguish  roughly  between  local  regions  where  a 
wave  with  a particular  frequency  oj  and  horizontal  wavelength  k , is  vertically  propa- 
gating ( ./'  > 0 ) or  vertically  evanescent  (/< 0).  We  shall  use  this  approach  for  our 
penumbral  model  in  Section  4. 

The  simplest  case  to  study  is  that  of  constant  c2,  t>i,  and  H.  We  shall  consider  this 
case  now  in  order  to  show  the  kinds  of  wave  modes  which  can  occur.  Although  this 
case  does  not  apply  strictly  to  a real  penumbra,  or  in  fact  to  our  penumbral  model,  we 
can  nevertheless  apply  the  resulting  dispersion  relation  locally  to  get  an  approximate 
picture  of  the  wave  modes.  We  can  also  approximate  a continuous  vertical  variation 
of  the  parameters  in  the  penumbra  by  a series  of  layers  in  which  they  are  constant.  In 
the  case  of  constant  c2,  v\.  and  H , Equations  (3)— (5)  have  constant  coefficients,  and 
we  can  assume  a solution  of  the  form  u(z)  = u exp(ik.z  + z/2H),  where  fi=  ( ii , f,  w)  is 
a constant  vector.  Here,  the  factor  exp(z/2W)  accounts  for  the  fact  that,  to  conserve 
energy,  the  perturbation  amplitude  must  grow  as  the  density  decreases.  The  system  of 
Equations  (3)-(5)  then  becomes  (cf..  Yu,  1965): 


(w2  - clk\)  u + ikx  |^/A.c2  - 

it-  = 0 , 

(12) 

(a>2  — v^k2x)  C = 0 , 

(13) 

r / c2  \i 

ikx  1 ikzc2  + ^</  - Jj  " + ~ ' 

/ , 1 ^ 
x ( (.  + . 

V ‘ 4 H2) 

) vv  — 0 . 

(14) 

The  dispersion  relation  for  waves  is  obtained  from  the  condition  for  nonzero  solu- 
tions of  the  homogeneous  system  ( 1 2)— ( 1 4),  i.e.,  the  vanishing  of  the  determinant  of 
the  coefficients.  This  yields  the  dispersion  relation 

(uj2  - v^k2)  jo/  - (c2  + l’a)  (aa2  + k2  + “2  + 

+ c2vik2x  + k*  + - 9 [9  ~ ^ j = 0-  (15) 

We  now  discuss  the  various  wave  modes  given  by  this  dispersion  relation,  with  some 
comments  on  their  relation  to  running  penumbral  waves. 
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2.1.  Thu  pure  acoustic  modi 

The  dispersion  relation  (15)  is  satisfied  with  u>2  = c2k2  for  the  particular  imaginary 
' alue  of  k.  given  by  ik.  = (\  2H)—  (ylc2).  From  Equations  ( !2)-(  14)  we  see  the  cor- 
responding motion  has  ;7#0,  whereas  f = w'  = 0.  This  mode  is  thus  a purely  compres- 
siona).  acoustic  mode  with  dispersion  relation  identical  to  that  of  a homogeneous 
isothermal  gas.  This  can  occur  since  the  motion  is  horizontal  and  parallel  to  the 
magnetic  field,  and  thus  there  is  no  contribution  from  buoyant  or  magnetic  forces. 
The  amplitude  behaves  as 


and  thus  grows  exponentially  with  height  (since  H <c2lg  for  stability  of  the  unpertur- 
bed atmosphere  (Yu.  1965)).  The  total  momentum  and  energy  are  finite,  however, 
provided  there  is  a lower  boundary  confining  the  motion  to  a semi-infinite  range  of  z. 
This  mode  is  identical  to  the  Lamb  mode  in  the  non-magnetic  case,  except  that  here 
the  scale  height  H is  modified  by  the  magnetic  field.  Since  this  mode  has  no  vertical 
motion,  it  can  not  be  associated  with  the  running  penumbral  waves. 

2.2.  The  pure  alfven  mode 

A root  of  the  dispersion  relation  (15),  for  arbitrary  k,,  is  given  by  u>2  = i >^k2,  which  is 
the  same  as  the  dispersion  relation  for  a pure  Alfven  wave  in  a homogeneous  atmos- 
phere with  uniform  magnetic  field.  With  w2  = v\k\,  however,  the  other  factor  ir.  the 
dispersion  relation  (15)  is  in  general  not  zero,  and  thus  Equations  (12)-(  14)  show  that 
the  motion  has  f^O.  whereas  u—  m =0.  The  motion  is  purely  horizontal  and  purely 
transverse  to  the  magnetic  field.  Thus,  the  pure  Alfven  mode  in  a stratified  atmosphere 
is,  so  to  speak,  plane  polarized,  with  no  vertical  motions.  With  vv=0  there  is  no  con- 
tribution from  the  buoyancy  force,  and,  further,  with  u= 0 the  motion  is  incompres- 
sible, so  the  wave  behaves  as  a pure  Alfven  wave  (with  the  amplitude  factor 
exp(z/2//)).  Again,  since  this  mode  has  no  vertical  motions,  it  can  not  be  associated 
with  the  running  penumbral  waves.  We  turn  now  to  the  remaining  roots  of  the  dis- 
persion relation  (15).  which  do  permit  vertical  motions. 

2.3.  The  plus  and  minus  modes 

The  remaining  roots  of  the  dispersion  relation  (15)  are  given  by 

»*.  - i {«’  + di  (d + d + Jp)  ± (V  + iff  (d  + d + 4;,y  - 
- ii-vid  [d  - d + Jy  - (o  - cHjdj  1.  06) 

We  shall  refer  to  these  two  modes  as  the  ‘plus’  and  ‘minus’  modes.  These  modes 
involve  the  interaction  of  all  three  restoring  forces:  buoyant,  pressure,  and  magnetic. 
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They  reduce  to  more  familiar  wave  modes  in  certain  limiting  cases.  For  example,  in 
the  limit  of  vanishing  magnetic  field  (t>A  -0)  they  reduce  to  the  acoustic  and  gravity 
modes  of  an  isothermal  atmosphere  (see,  for  example,  Thomas  cr  at.,  1971).  Alter- 
natively, in  the  limit  of  no  stratification  (//  —•  x,  y ->0),  they  reduce  to  the  fast  and 
slow  magnetoacoustic  waves  in  a homogeneous  atmosphere  (see,  for  example,  Oster- 
brock,  1961 ).  For  intermediate  cases  (such  as  for  penumbral  conditions),  we  can  look 
at  the  modes  as  being  either  magnetically  modified  acoustic-gravity  waves  or  gravity- 
modit.ed  magnetoacoustic  waves.  We  shall  continue  to  use  the  terms  plus  and  minus 
modes  here. 

The  plus  and  minus  modes,  for  k,=  0,  are  shown  schematically  in  a diagnostic 
diagram  in  Figure  1.  The  asymptotic  behavior  of  the  dispersion  relation  (16)  for 
rrt  = 0 is  such  that 


, c + cA  , 

UJ+  -*•  ...i  > w-  - 0 lor  kx  ->  0, 

4/7 

and 

toi  - k2x  max(c2,  e2),  co2_  - k2x  min  (c2,  t>A)  for  kx->oo. 


(17) 

(18) 


Fig.  1.  Schematic  diagnostic  diagram  of  the  plus  and  minus  modes,  showing  curves  for  k,  ^ 0 
(solid  lines).  The  plus  mode  has  a finite  cutoff  frequency  as  £*->-0.  For  kx-+cc,  the  plus  and  minus 
modes  approach  the  dispersion  lines  for  the  pure  acoustic  and  pure  Atfven  modes.  For  strong  stratifi- 
cation and  weak  magnetic  field,  the  minus  mode  approaches  the  line  c»~kxi\  from  above  rather 

than  below. 
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lu  the  diagnostic  diagram  (figure  I).  there  are  vertically-propagating  waves  (A;  >0) 
above  the  m,  curve  and  below  the  <u  curve,  and  evanescent  waves  (A2<0)  in  the 
region  between  the  two  curves. 

The  plus  mode  has  a finite  cutoff  frequency  of  (wt),=  (c2  + vji)',2l2H  as  A\,-*0. 
Estimates  of  this  cutoff  frequency  at  photospheric  levels  in  a penumbra  give  values 
just  a little  lower  than  the  frequencies  of  penumbral  waves.  This  suggests  that  the 
penumbral  waves  might  be  identified  with  plus  modes.  We  shall  show  this  in  more 
detail  through  the  use  of  our  penumbral  model. 


3.  The  Penumbral  Model 


We  now  present  the  penumbral  model  to  be  used  in  our  wave  calculations.  For  our 
purposes,  a penumbral  model  consists  of  specified  distributions  of  c2,  v£,  and  H with 
height  z.  To  simplify  the  calculations,  we  have  chosen  to  represent  the  expected  verti- 
cal variation  of  these  parameters  by  piecewise  linear  functions.  The  specified  forms 
of  c2  (z),  i>a  (- ),  and  H (- ) are  shown  in  Figure  2.  The  important  features  of  the  model 
are  as  follows.  The  sound  speed  increases  with  depth  into  the  convection  zone  from 
a broad  minimum  in  the  penumbral  photosphere  and  chromosphere.  The  Alfven 
velocity  increases  rapidly  with  height  due  to  the  nearly  exponential  decrease  of 
density  while  the  magnetic  field  decreases  nearly  linearly.  The  density  scale  height 
increases  on  either  side  of  a minimum  in  the  penumbral  photosphere  and  low  chromo- 
sphere. The  model  thus  reproduces  the  main  expected  features  of  vertical  penumbral 
structure.  The  numerical  values  of  the  parameters  were  chosen  to  represent  a typical 
penumbra. 

The  distribution  of  the  sound  speed  was  determined  primarily  from  the  expected 
temperature  distribution.  Kjeldseth  Moe  and  Maltby  (1969)  report  that  the  tempera- 
ture T in  the  penumbra  may  be  obtained  by  adding  a constant  <40  = 0.055  to  the  0 
values  of  the  quiet  photosphere,  where  0 = 5040/7.  With  this  AO,  the  calculated  relative 
intensities  averaged  over  the  penumbral  fine  structure  agreed  well  with  observations. 
Kjeldseth  Moe  and  Maltby  used  the  Bilderberg  Continuum  Atmosphere  (Gingerich 
and  de  Jager,  1968)  for  their  quiet  photospheric  temperatures;  here,  we  use  the  more 
recent  Harvard-Smithsonian  Reference  Atmosphere  (HSRA,  Gingerich  el  at.,  1971). 
Using  the  constant  AO =0.055,  the  penumbral  temperature  minimum  is  found  to  be 
3989 K,  yielding  a minimum  value  of  the  sound  speed  squared  of  c2  = 43.5  km2  s'2. 
The  slope  of  the  linear  increase  of  c2  into  the  qonvection  zone  was  taken  to  be  0.1 
km  s'2;  this  choice  is  also  based  on  the  behavior  of  the  HSRA.  The  distribution  of 
sound  speed  (in  km  s'1)  is  given  by 


c2(z) 


63.5  - O.lz 

43.5 


z < 200 
z > 200. 


(19) 


At  some  point  in  the  penumbral  chromosphere  the  sound  speed  will  increase  rapidly 
to  coronal  values.  However,  this  increase  takes  place  above  the  region  where  the 
penumbral  waves  are  trapped  in  our  model,  so  we  have  not  included  this  in  the  model. 
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2.10"*  (km) 


Fig  2.  Distribution  of  c2,  t A2,  and  H with  height  z in  the  pcnumbral  model  (solid  lines).  The  data 
(crotses)  for  rA2  are  based  on  penumbral  observations  (sec  text).  Data  points  for  quiet  photosphcric 
values  of  c2  and  H,  based  on  the  HSR  A.  arc  shown  only  for  comparison  with  the  penumbral  mocel. 

v 

The  distribution  of  cz  is  shown  in  Figure  2.  along  with  quiet  atmospheric  values  of  c1 
from  the  HSRA  (Nakagawa,  1973)  for  comparison. 

To  determine  the  Alfven  velocity,  we  assumed  that  the  magnetic  field  strength 
decreases  linearly  with  height  and  is  1000  G at  2 = 0.  The  rate  of  decrease  of  field 
strength  with  height  was  taken  to  be  0.2  G km'1,  in  accordance  with  observations  (see 
Bray  and  Loughhead,  1964).  The  photosphcric  densities  were  estimated  from  the 
penumbral  model  of  Makita  (1963),  the  only  penumbral  model  to  give  densities.  His 
densities  are  consistent  with  more  recent  quiet  photospheric  models.  The  resulting 
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values  of  t*  are  plotted  in  Figure  2.  The  \ertical  distribution  of  these  points  was 
approximated  in  our  model  by  two  straight-line  segments,  with  the  upper  seemcnt 
extrapolated  to  greater  heights.  In  the  convection  zone,  the  Alfven  velocity  was 
approximated  by  another  straight  line  segment  which  goes  to  zero  at  ; = -500  km, 
due  primarily  to  the  increase  in  density  with  depth,  but  also  to  the  fact  that  the  penum- 
fcral  field  probably  lies  over  the  convection  zone.  The  resulting  distribution  of  the 
Alfven  velocity  (in  km  s ')  is  given  by 


i A (2) 


0 

25  -f-  0.05; 

25  +0.172; 

94  + 0.75  (z  - 400) 


; < - 500 
- 500  « z < 0 
0 sj  z « 400 
z > 400. 


(20) 


The  distribution  of  the  local  density  scale  height  H (in  km)  with  height  was  chosen 
to  be 

I 230  - 0.375;  ; < 200 

//(;)=  155  200  ^z  <400  (21) 

I 155  +0.55  (z  - 400)  z^400. 

In  choosing  values  for  //,  we  were  guided  by  the  following  expression  for  the  vertical 
ertropy  gradient  in  the  atmosphere: 


da 

OC,K 

1 1 d 1 

f"*V| 

d; 

= PT 

LV  2; 

1 H d;  1 

U/J 

where  r,.  is  the  specitic  heat  at  constant  volume,  k is  the  isothermal  compressibility, 
and  /f  is  the  coefficient  of  thermal  expansion.  This  relation  follows  from  the  hydro- 
static equation  1 1 ) and  the  basic  thermodynamic  relation 


Tds  = 


2 

(dp  — c dp). 


(23) 


The  distribution  of//  given  in  (21 ) is  such  that  the  entropy  gradient  d.v/dz  is  negative  in 
the  upper  convection  zone  and  positive  in  the  penumbral  photosphere  and  chromo- 
sphere. The  distribution  of  //  is  shown  in  Figure  2,  along  with  quiet  atmospheric 
values  of  //  from  the  FISRA  (Nakagawa.  1973)  for  comparison.  Flere  we  can  see  the 
effect  of  the  supporting  magnetic  field;  the  scale  height  is  greater  in  the  penumbra 
than  in  the  quiet  atmosphere.  The  effect  of  the  magnetic  field  increases  as  the  density 
decreases. 

4.  Trapping  of  Plus  Modes  in  the  Penumbral  Photosphere 

We  now  show  that  the  running  penumbral  waves  can  be  identified  in  our  model  with 
plus  modes  that  are  vertically  trapped  at  photospheric  levels.  The  basic  mechanism  for 
the  trapping  is  the  refraction  due  to  (i)  the  increasing  Alfvin  velocity  with  increasing 
height  in  the  photosphere-low  chromosphere,  and  (ii)  the  increasing  sound  speed  with 
deptn  in  the  convection  zone.  We  shall  demonstrate  the  trapping  mechanism  in  two 
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ways:  lirst,  by  considering  the  local  dispersion  relation  (16)  at  three  distinct  levels  in 
the  model  penumbra  (assuming  locally  constant  parameters),  and  second,  by  con- 
sidering the  propagation  equation  (10)  for  the  complete  continuous  pcnumbral 
model  (Figure  2). 

4.1.  Three-level  model 

A convenient  method  of  looking  at  wave  modes  in  our  model  is  to  draw  diagnostic 
diagrams  such  as  Figure  1 for  various  heights  in  the  atmosphere,  assuming  locally 
constant  values  of  c2,  vA,  and  H - that  is,  ignoring  the  derivatives  of  these  parameters 
in  Equation  (3)-(5).  Although  only  approximate,  this  method  does  give  some  feeling 
for  the  behavior  of  the  wave  modes.  We  shall  draw  diagnostic  diagrams  for  three 
different  levels  in  our  model  penumbra;  one  in  the  convection  zone,  one  in  the 
photosphere,  and  one  in  the  chromosphere.  We  may  interpret  these  diagnostic 

\«10"5(krn) 


k*  »10J  (km"1) 

Fig.  3.  Superimposed  diagnostic  diagrams  for  values  of  c1,  vA2,  and  H at  three  different  levels  in 
the  model  penumbra : level,  I r=  - 1400  km;  level  2,  : — — 50  km;  level  3,  r - 500  km.  The  shaded 
region  is  a region  of  vertical  trapping  of  plus  modes  around  level  2.  The  crosses  correspond  to 
observed  running  penumbral  waves  (Giovanelli,  1974). 
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diagram;,  more  properly  as  representing  a three-layer  model  of  the  penumbra,  where 
in  each  layer  the  parameters  are  constant.  A middle  layer  of  Unite  thickness  represent- 
ing the  pcnumbral  photosphere  is  bounded  above  and  below  by  semi-infinite  layers 
representing  the  penumbral  chromosphere  and  convection  zone.  respectively. 

Figure  3 shows  diagnostic  diagrams  for  values  of  c2.  r*.  and  //  at  three  different 
levels  in  the  penumbral  model:  r = - 1400  km  (layer  I,  convection  /one),  z = -50  km 
(layer  2.  penumbral  photosphere),  and  ’ = 500  km  (layer  3,  penumbral  chromosphere). 
The  important  feature  of  this  figure  is  the  existence  of  the  shaded  region  in  which  the 
plus  modes  are  vertically  propagating  in  layer  2,  but  are  vertically  evanescent  in 
layers  I and  3.  These  modes  are  thus  trapped  in  layer  2.  The  plus  mode  diagnostic 
curve  for  layer  3 lies  above  that  of  layer  2 because  of  the  higher  Alfven  velocity  in 
layer  3.  For  higher  levels  in  the  chromosphere,  the  plus  mode  curve  for  layer  3 will  be 
higher  than  that  shown.  The  plus  mode  diagnostic  curve  for  layer  I has  a lower  cutoff 
frequency  than  layer  2 due  to  the  larger  scale  height  in  the  convection  zone,  but  lies 
above  the  curve  for  layer  2 for  higher  values  of  kx  because  of  the  higher  value  of  c2  in 
layer  I . For  deeper  levels  in  the  convection  zone,  the  plus  mode  curve  for  layer  1 will 
have  an  even  lower  cutoff  frequency  and  a steeper  slope.  The  existence  of  the  region  of 
trapping  in  the  diagnostic  diagram  is  a consequence  of  the  qualitative  features  of  our 
model  and  is  not  dependent  on  the  particular  choice  of  numerical  values. 

There  is  no  trapping  of  minus  modes  in  Figure  3.  All  the  minus  modes  propagating 
in  layer  2 are  also  propagating  in  layer  3.  The  minus  mode  curve  for  layer  ! has 
vanished  since  the  Alfven  velocity  is  zero  and  the  atmosphere  is  convectively  unstable 
at  that  level. 

Also  shown  in  Figure  3 are  four  data  points  corresponding  to  penumbral  waves 
observed  in  different  sunspots  by  Giovanelli  (1974)  for  which  he  gives  specific  wave- 
length and  periods.  These  points  tend  to  cluster  in  the  long-wavelength  end  of  the 
shaded  region  of  trapping  of  the  plus  modes.  The  one  point  which  lies  outside  the 
shaded  region  does  lie  in  the  region  of  trapping  if  the  parameter  values  in  layer  I are 
chosen  to  represent  a lower  level  in  the  convection  zone.  That  is,  a wave  of  this  fre- 
quency and  wavelength  is  reflected  at  a lower  level  in  the  convection  zone  in  our 
model.  This  data  point  is  considered  atypical  by  Giovanelli,  however.  The  observation- 
al data  correspond  to  different  sunspots,  no  doubt  having  different  field  strengths, 
wnereas  the  diagnostic  curves  are  for  a single  choice  of  the  model  parameters.  Never- 
theless, the  resulting  picture  in  Figure  3 clearly  shows  that  the  running  penumbral 
waves  should  be  identified  with  plus  modes  which  are  vertically  trapped  at  photo- 
spheric  levels  in  our  model. 

4.2.  Continuous  model 

We  now  illustrate  more  accurately  the  trapping  of  the  penumbral  waves  by  making  use 
of  our  complete  penumbral  model  (Figure  2)  in  conjunction  with  the  propagation 
equation  (10).  Roughly  speaking,  a wave  is  vertically  propagating  when  the  function 
/ tu;  z)  (Equation  (1 1))  is  positive,  and  is  vertically  evanescent  when/is  negative. 
In  Figure  4 we  have  plotted /as  a function  of  z for  our  penumbral  model  (Figure  2)  for 
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Fig.  4.  The  function  / ( kx , to;  ;)  evaluated  for  the  pcnumbral  model  (Figure  2)  with  kx  - 2 x I0~ 3 

km'1,  oj  = 3 x 10"  2 s'1. 

a typical  value  of  horizontal  wavenumber  (kx  = 2 x 10~3  km"1 ) and  frequency 
(w  = 3x  I0~2  s"1)  for  penumbral  waves.  The  function  / (z)  is  discontinuous  due  to 
the  discontinuities  in  the  first  derivatives  of  the  parameters  c2,  and  H in  our  model. 
The  general  behavior  of  / is  that  it  is  positive  in  a central  region  extending  from 
z*  — 1000  km  to  roughly  z»300  km,  and  negative  above  and  below  this  region. 

In  the  lowest  order  WKB  approximation,  the  general  solution  to  (10)  is  given  by 


<t>(z)  = exp  / J J f dzj  . (24) 

o 

If /is  positive  for0<z<z0  and  negative  for  z>z0,  then  the  solution  for  z>z0  which  is 
bounded  as  z-*co  can  be  written  as 

Zo  z 

<p(z)  = exp  ± i j"  Jfdz- j J -/dzj.  (25) 

O 20 

Thu>,  for  z>z0,  the  solution  has  the  form  of  an  exponential  decay  (with  variable 
exponential  factor),  and  the  wave  may  be  characterized  as  evanescent  in  the  vertical 
direction.  This  analysis  is  somewhat  crude,  however,  since  the  WKB  solution  (24)  is 
not  accurate  in  the  neighborhood  of  z0.  The  same  analysis  applies  to  the  case  where/ 
is  positive  for  -z|,<z<0  and  negative  for  z<  — z'0.  The  wave  then  becomes  evanes- 
cent as  z decreases  below  -z'0.  Applying  this  to  Figure  4,  we  conclude  that  a wave 
with  the  specified  horizontal  wavelength  and  frequency  is  trapped  roughly  in  the 
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central  region  — I000<z<300  where /is  positive.  The  exact  location  of  the  upper 
boundary  of  the  trapping  region  is  somewhat  unclear  due  to  the  complicated  dis- 
continuous behavior  off  there.  This  is  due  to  the  piecewise  linear  nature  of  our  model; 
with  smoother  distributions  of  c1 , v\,  and  H , the  function  / would  have  a single, 
smooth  zero  crossing  at  some  z>(J. 

We  have  been  discussing  the  behavior  of  the  function  4>{z).  Comparison  with  the 
case  of  constant  parameters  shows  that  |0j2  is  roughly  proportional  to  the  kinetic 
energy  of  vertical  motion  of  the  wave  w).  Thus,  it  is  the  kinetic  energy  that  is 

trapped  in  the  region  of  positive  /.  The  amplitude  of  the  vertical  velocity  behaves 
somewhat  differer.  From  Equations  (9)  and  (25)  we  have,  for  z>z0. 


«•(-)=  <•  'i±l  j v//dz  + C(2) 


( - ) = ~ i J A dz  - J J - f dr . 


(26) 

(27) 


The  function  A (z)  (Equation  (7))  is  negative  over  a range  in  z,  and  the  function  Q(z) 
will  remain  positive  over  a range  of  z beyond  z0  before  becoming  negative.  Thus,  the 
vertical  velocity  continues  to  increase  above  z0  before  reaching  a maximum  at  z = z,, 
where  z,  is  determined  from  the  relation  Q(z ,)  = 0.  This  can  occur  in  spite  of  the 
decreasing  energy  above  z = zu  due  to  the  rapid  decrease  in  density. 

By  numerically  evaluating  the  integrals  in  (27),  we  have  estimated  the  point  z,  of 
maximum  vertical  velocity  in  our  model,  again  for  the  choice  kx  = 2 x 10' 3 km'1  and 
tu  = 3 x )0'2  s"1.  The  resulting  value  is  z,  « 1250  km.  This  is  at  a level  where  the  wave 
is  evanescent,  and  is  well  into  the  region  of  formation  of  Ha  (Vernazza  et  al.,  1973). 
Thus,  although  the  trapped  waves  have  their  maximum  energy  at  lower  levels  (z^O) 
where  the  density  is  higher,  the  vertical  velocity  is  greatest  at  chromospheric  levels 
where  is  it  observed. 


5.  Discussion 


An  important  conclusion  that  emerges  from  our  model  is  that  the  running  penumbral 
waves  are  basically  a photospheric  phenomenon,  even  though  they  are  observed  at 
chromospheric  levels.  This  is  really  an  expected  result;  the  phase  velocity  of  the 
penumbral  waves  is  typical  of  photospheric  conditions  rather  than  chromospheric. 
For  example,  the  Alfven  velocity  at  Ha  levels  in  the  penumbra  is  far  greater  than  the 
phase  velocity  of  the  penumbral  waves.  In  our  picture,  the  waves  are  evanescent  where 
they  are  observed.  That  is,  the  observed  vertical  motions  are  a passive  response  to  an 
actively  propagating  wave  at  lower  levels.  We  have  shown  that  although  the  energy  is 
trapped  at  photospheric  levels,  the  vertical  velocity  reaches  its  maximum  value  at 
levels  of  formation  of  Ha,  consistent  with  observations. 

The  concept  of  the  penumbral  waves  as  a photospheric  phenomenon  is  also  quite 
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consistent  with  their  expected  source  of  excitation.  Moore  (1973)  has  shown  that  the 
likely  source  of  excitation  of  the  penumbral  waxes  is  oscillatory  convection  in  a sub- 
photospheric  layer  in  the  umbra.  The  observation  that  the  penumbral  waves  dis- 
appear suddenly  at  the  boundary  between  the  white-light  penumbra  and  the  sur- 
rounding photosphere  also  supports  the  photospheric  nature  of  the  waves.  In  the 
chromosphere,  the  penumbral  fibril  structure  extends  outward  beyond  this  boundary. 

The  basic  mechanism  for  the  vertical  trapping  of  the  penumbral  waves  is  not  sensi- 
tive to  the  details  of  our  penumbral  model.  I his  mechanism  is  based  primarily  on  the 
increase  in  the  sound  speed  as  we  go  down  into  the  convection  /.one  and  the  increase 
in  the  Alfven  velocity  as  we  go  up  into  the  chromosphere. These  features  are  certain 
to  remain  in  any  improved  penumbral  model.  Observationally.  the  insensitivity  of  the 
basic  mechanism  is  confirmed  by  the  fact  that  penumbral  waves  are  seen  in  almost 
every  stable,  regular  sunspot  (Zirin  and  Stein.  1972). 

The  character  of  the  trapped  plus  modes  is  different  at  different  heights.  The  wave  is 
more  nearly  acoustic  at  low  levels  (convection  zone),  but  is  more  nearly  Alfvenic  at 
higher  levels  (photosphere  and  low  chromosphere).  The  effect  of  stratification  and 
gravity  , while  not  dorr  inant.  is  felt  throughout  the  trapping  region. 

Although  vve  have  not  attempted  to  calculate  them,  there  are  certainly  resonant 
modes  in  our  model  which  arise  from  the  constructive  interference  of  plus  modes 
reflected  from  above  and  below  in  the  trapping  region.  The  resonant  dispersion  curves 
would  lie  in  the  shaded  region  in  Figure  3.  However,  since  the  penumbral  waves  are  a 
somewhat  transitory  phenomenon,  there  is  no  reason  to  associate  them  specifically 
with  such  a resonant  mode:  it  is  sufficient  that  a wave  of  the  proper  frequency  and 
horizontal  wavelength  be  vertically  trapped. 

Since  our  penumbral  model  is  horizontally  uniform,  it  cannot  account  for  any 
observed  horizontal  (radial)  variations  in  the  penumbral  waves.  We  can  roughly 
account  for  the  radial  geometry  of  the  penumbra  by  noting  that  the  energy  in  a verti- 
cally trapped  wave  will  decrease  as  I jr,  where  r is  radius  measured  from  the  spot 
center,  as  the  wave  propagates  radially  outward.  This  accounts  at  least  in  part  for  the 
observed  decrease  in  wave  amplitude  with  radius.  The  sudden  disappearance  of  the 
waves  at  the  outer  edge  of  the  penumbra  is  associated  with  the  more  fundamental 
problem  of  the  basic  existence  of  the  sharp  boundary  between  the  white-light  penum- 
bra and  the  surrounding  photosphere. 

Finally,  we  should  note  that  our  assumption  of  adiabatic  perturbations  is  invalid 
over  a limited  range  of  height  in  the  low  photosphere,  say  the  first  few  hundred 
kilometers  above  c = 0.  Here  one  should  account  for  the  rapid  radiative  exchange.  The 
main  effect  of  the  radiative  transfer  on  the  plus  modes  is  to  effectively  replace  the 
adiabatic  sound  speed  in  this  limited  region  by  something  nearer  to  the  isothermal 
sound  speed.  The  overall  effect  on  the  trapped  plus  modes  discussed  here  w ill  be  quite 
small,  even  quantitatively. 
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Convective  instability  in  the  presence  of  a nonuniform 
horizontal  magnetic  field 
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Newcomb's  criterion  for  convective  stability  in  the  presence  of  a horizontal  magnetic  field  is  written 
in  a form  which  explicitly  shows  the  effect  of  vertical  variations  of  the  magnetic  field  strength.  It  is 
shown  that  a nonuniform  horizontal  magnetic  field  can  be  destabilizing  as  well  as  stabilizing 


We  consider  the  convective  instability  of  a compres- 
sible, inviscid,  perfectly  conducting  gas  permeated  by 
a horizontal  magnetic  field  B = [/?(*),  0,  0]  that  may  vary 
with  height  z,  under  a uniform  gravitational  accelera- 
tion g (in  the  negative  z direction),  assuming  adiabatic 
perturbations.  The  first  complete  treatment  of  this 
problem  was  given  by  Newcomb1  using  the  energy  inte- 
gral method.  He  showed  that  a necessary  and  sufficient 
condition  for  stability  is  given  by 


dp  >pV 

dz  yp  ’ 


(1) 


where  p is  mass  density,  p is  pressure,  and  y is  the 
ratio  of  specific  heats.  In  the  case  of  instability,  the 
most  unstable  mode  has  the  form  of  an  interchange  of 
long  but  finite  segments  of  magnetic  field  lines.  The 
stability  criterion  (1)  was  also  derived  by  Yu2  by  con- 
sidering the  force  balance  on  a displaced  magnetic  flux 
tube. 

Newcomb  noted  that  the  critical  density  gradient  on 
the  right-hand  side  of  (1)  is,  at  least  explicitly,  inde- 
pendent of  the  magnetic  field,  and  is,  in  fact,  that  given 
by  the  Schwa rzschild  criterion3  in  the  absence  of  a mag- 
netic field.  However,  in  the  case  of  a nonuniform  hori- 
zontal magnetic  field  (for  which  Newcomb’s  analysis  is 
valid),  the  static  distribution  of  pressure  and  density  is 
affected  by  the  magnetic  field,  and  thus  the  stability 
criterion  (1)  depends  implicitly  on  the  magnetic  field. 
This  point,  which  was  overlooked  by  Newcomb,  will  be 
pursued  here. 


Using  the  equation  of  hydrostatic  equilibrium, 


and  the  equation  of  state  of  a perfect  gas,  p = pRT,  the 
stability  criterion  (1)  can  be  written  in  the  form 


dT  _(dT\  > J_  d/fifi 
dz  \dz  /,  pR  dz  \8rr/ 


(3) 


where  (dT/dz)s  = -g/c,  = -g(v-  \)/yR  is  the  adiabatic 
temperature  gradient.  The  new  form  (3)  of  the  stability 
criterion  has  an  advantage  over  the  form  (1)  in  that  it 
shows  explicitly  the  effect  of  vertical  variations  of  the 
magnetic  field,  since  the  temperature  distribution  is 
independent  of  the  magnetic  field.  For  a uniform  mag- 
netic field  (J5=  const),  (3)  reduces  to 


which  is  identical  to  the  Schwarzschild  criterion  in  the 
absence  of  a magnetic  field.  [Note  that  usually,  e.g. , 
in  a star,  dT/dz  and  ( dT/dz )s  are  both  negative.]  Thus, 
a uniform  horizontal  magnetic  field  has  no  effect  on  the 
condition  for  the  onset  of  convective  instability, 4 al- 
though, as  Newcomb1  also  showed,  it  does  have  an  ef- 
fect on  the  growth  rates  of  unstable  modes. 


If,  however,  the  magnetic  field  is  nonuniform,  then 
(3)  shows  that  the  field  can  be  stabilizing  (in  the  case 
rlB/dz  >0)  or  destabilizing  ( dB/dz  <0).  A field  that  in- 
creases with  height  ( dB/dz  >0)  can  stabilize  the  atmo- 
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sphere  in  the  presence  of  a superadiabatlc  temperature 

gradient.  When  .he  field  decreases  with  height  ( tIB/dz 
sO),  the  critical  temperature  gradient  is  reduced  be- 
low the  adiabatic  gradient. 


As  an  illustration  of  the  destabilizing  effect  for  dB/ 
dz  <0,  consider  the  case  of  an  isothermal  at  nosphere, 
which  is  convectively  stable  in  the  absence  of  a mag- 
netic field.  The  stability  criterion  (3)  reduces  to 


If  the  magnetic  piessure  decreases  more  raridly  with 
height  than  the  critical  rate 


the  atmosphere  will  be  convectively  unstable.  For  this 
critical  magnetic  field  gradient,  the  corresponding 
pressure  gradient  is  found  from  (2)  to  be 


Thus,  in  the  state  of  marginal  stability  in  an  isothermal 
atmosphere,  the  gradient  of  gas  pressure  balances  the 
fraction  l/y  of  the  gravitational  force,  while  the  grad- 
ient of  magnetic  pressure  balances  the  remaining  frac- 
tion (y  - 1 )/y  of  the  gravitational  force. 

A special  case  of  interest  is  that  of  an  isothermal  at- 
mosphere in  which  the  field  decreases  with  height  in 
just  the  manner  that  causes  the  Alfven  speed  i>A  = (B2/ 
4irp)1/2  to  remain  constant.  In  this  case,  (4)  can  be 
written  as 


where  r = (yRT)xn  is  the  adiabatic  sound  speed.  This 
stability  criterion  was  obtained  for  this  special  case  by 
Yu5  from  a normal  mode  analysis.  Parker®  has  studied 
this  case,  including  the  effect  of  cosmic -ray  pressure, 
in  connection  with  the  gaseous  disk  of  the  galaxy. 

The  destabilizing  effect  of  a nonuniform  magnetic 
field  with  dB/dz  <0  is  similar  to  the  phenomenon  of 
“magnetic  buoyancy”  analyzed  by  Parker7  and  by  Weiss.® 
Magnetic  buoyancy  is  attributed  to  an  isolated  tube  of 
magnetic  flux  in  thermal  equilibrium  with  its  nonmag- 
netic surroundings.  The  instability  discussed  here  is 
due  to  the  buoyancy  of  an  arbitrary  tube  of  flux  within 
a smoothly  varying  magnetic  field. 
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ABSTRACT 

The  linearized  theory  of  magneto-atmospheric  waves  (involving  the  combined  restoring  forces 
due  -o  buoyancy,  compressibility,  and  magnetic  iic.d)  is  developed  for  the  case  of  a horizontal 
magnetic  held.  A general  propagation  equation  is  derived  for  adiabatic  perturbations  with  arbitrary 
vertical  distributions  of  the  sound  speed  r,  Alfven  velocity  i\,  and  local  density  scale  height  H. 

An  exact  analytical  solution  to  the  propagation  equation  is  obtained  for  the  case  of  an  isothermal 
atmosphere  permeated  by  a uniform  horizontal  magnetic  held,  without  making  the  usual  short- 
wavelength  assumption.  This  solution  is  applied  to  an  idealized  model  of  the  low-corona- 
chromosphere  transition  region  for  comparison  w ith  observations  of  flare-induced  coronal  waves. 

The  results  show  that  disturbances  may  propagate  horizontally  in  the  low  corona  in  a wave 
guide  formed  by  the  sudden  density  increase  into  the  chromosphere  below  and  by  the  rapidly 
increasing  Alfven  velocity  with  height  in  the  corona.  The  group  velocity  of  the  guided  wave  modes 
is  nearly  independent  of  wavelength,  so  that  a disturbance  propagates  as  a compact  wave  packet. 

Subject  headings:  hydromagnetics  — Sun:  atmospheric  motions  — Sun:  corona  — • 

Sun:  magnetic  fields 

I.  INTRODUCTION 

The  theory  of  waves  in  a compressible,  stratified,  electrically  conducting  atmosphere  permeated  by  a mag- 
netic field  is  of  considerable  importance  in  astrophysics,  especially  in  solar  physics  where  there  is  a wealth  of  detailed 
observations  of  such  waves  in  the  solar  atmosphere.  Following  Yu  (1965),  we  shall  refer  to  waves  which  involve 
compressibility,  buoyancy,  and  magnetic  forces  as  magneto-atmospheric  naves.  Among  the  many  solar  phenomena 
that  are  seemingly  attributable  to  magneto-atmospheric  waves  are  the  heating  of  the  corona  and  of  chromospheric 
plages,  the  5-minute  oscillations  in  active  regions,  oscillations  in  sunspot  umbrae  and  penumbrae,  and  flare- 
induced  coronal  disturbances. 

The  theory  of  magneto-atmospheric  waves  is  complicated  by  the  anisotropic  nature  of  the  medium;  the  gravita- 
tional field  and  the  magnetic  field  each  introduce  preferred  directions.  Additionally,  the  disturbance  is  subjected 
to  the  combined  restoring  forces  due  to  compressibility,  buoyancy,  and  the  magnetic  field,  so  that  pure  wave 
modes  (i.e.,  acoustic,  Alfven,  and  gravity)  exist  only  as  special  cases.  In  general,  magneto-atmospheric  waves 
involve  the  effects  of  all  three  restoring  forces.  The  problem  is  compounded  by  the  fact  that  the  basic  parameters 
describing  wa«  propagation  in  the  solar  atmosphere  (the  sound  speed  c,  the  Alfven  velocity  t A,  and  the  local 
density  scale  height  H)  are.  in  general,  functions  of  height,  and  therefore  the  disturbance  cannot  be  represented  by 
plane  waves  propagating  in  the  atmosphere. 

Previously,  the  problem  of  magneto-atmospheric  waves  has  been  studied  by  either  of  two  basic  approaches,  each 
yielding  a dispersion  relation  based  on  constant  values  of  the  atmospheric  parameters.  The  first  approach  has  been 
to  assume  that  the  vertical  extent  of  the  disturbance  is  much  less  than  the  smallest  scale  height  for  variation  of  the 
atmospheric  parameters.  The  parameters  can  then  be  taken  as  constant  locally.  McLellan  and  Winterberg  (1968) 
studied  an  isothermal  atmosphere  permeated  by  a uniform  magnetic  field  with  arbitrary  orientation.  Then,  assum- 
ing that  the  Alfven  velocity  is  constant  locally  (although  it  actually  increases  exponentially  with  height),  they 
derived  a local  dispersion  relation  that  is  valid  for  short  wavelengths  (short  compared  with  the  density  scale  height). 
This  local  dispersion  relation  has  been  studied  by  several  authors  (Bel  and  Mein  1971;  Michalitsanos  1973; 
Nakagawa  et  ai.  1973)  to  determine^the  effects  of  different  propagation  directions  and  magnetic  field  orientations. 

The  second  basic  approach  has  been  to  investigate  an  isothermal  atmosphere  permeated  by  a horizontal  magnetic 
field  that  decreases  exponentially  with  height  such  that  the  Alfven  velocity  is  constant,  as  are  the  sound  speed 
and  the  density  scale  height.  Yu  (1965)  derived  the  dispersion  relation  whicji  in  this  case  is  valid  for  all  wavelengths. 
He  evaluated  the  three  modes  of  propagation  for  various  angles  between  the  wave  propagation  vector  and  the 
magnetic  field.  Chen  and  Lykoudis  (1972)  used  the  dispersion  relation  of  Yu  to  study  the  5-minute  oscillations  in 
plage  regions.  Nye  and  Thomas  (1974a)  used  Yu’s  dispersion  relation  in  connection  with  a multilayer  model  for 
running  penum  rral  waves. 

• Also  C.  E.  Kenneth  Mees  Observatory. 
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We  consider  only  plane-parallel  atmospheres,  with  no  horizontal  variation,  in  order  to  permit  Fourier  dc- 
coinposition  in  the  hori/ontal  directions  as  well  as  lime.  With  variations  in  the  .’-direction  only,  the  governing 
partial  dillerential  equations  reduce  lo  ordmarv  dilVcreniial  equations.  There  are  only  iwo  magnelic  liefd  con- 
figurations eoiisislcnl  with  sialic  equilibrium  ,unl  no  hori/ontal  variation.  I hey  arc  B constant  and  It 

/<„(;),  0).  We  shall  restrict  our  study  lo  the  case  of  a unidirectional  hori/ontal  magnetic  held  that  may  vary 
with  height.  We  shall  not  follow  either  of  I lie  Iwo  basic  approaches  discussed  above,  however,  because  of  the 
inherent  limitations  ol  r.u  h of  them  lor  solar  applications  I lie  length  scale  for  magnetic  field  changes  in  the  solai 
utmnsphrir  is  gcuei.dty  linn  li  gic.itn  Itiuu  tin  dricaly  stale  liriglil,  so  ilml  the  assumption  til  constant  Allvcli 
velocity  is  not  jt  stilted  < m the  other  hand,  the  length  scale  ol  observed  disturbances  in  ihe  solar  atmosphere  is  not 
generally  small  compared  with  the  density  scale  height.  This  is  especially  true  in  the  photosphere  where  the  density 
scale  height  may  be  smaller  than  100  km,  and  hence  smaller  than  the  limit  of  observational  resolution. 

In  § If  we  derive  a general  propagation  equation  for  an  arbitrary  direction  of  propagation  and  arbitrary  vertical 
variations  of  the  atmospheric  parameters.  For  the  case  of  an  isothermal  atmosphere  permeated  by  a uniform 
horizontal  magnetic  field,  the  Alfven  velocity  increases  exponentially  with  height.  In  §111  we  obtain  an  exact 
general  solution  of  the  propagation  equation  in  this  case.  We  compute  eigenmodes  for  the  case  of  a rigid  lower 
boundary  in  § IV,  and  apply  this  to  a specific  solar  wave  phenomenon,  the  flare-induced  coronal  waves,  in  § V. 

The  analytical  treatment  in  §§  11  and  111  also  forms  the  basis  of  a following  paper  in  which  we  deal  with  running 
penumbral  waves. 


II.  BASIC  EQUATIONS 

The  atmosphere  is  assumed  to  be  a compressible,  inviscid,  perfectly  conducting  gas  under  a uniform  acceleration 
of  gravity  # ( = 0 274  kmr2)  in  the  negative  r-direction.  The  undisturbed  magnetic  field  is  taken  in  the  x-direction 
and  may  vary  with  height  r;  i.c.,  B0  = [fi0(-).  0-  0].  The  undisturbed  pressure,  density,  and  temperature  may  all  be 
functions  of  height  r,  and  are  denoted  by  p0(z),  p0(z).  and  7~0(- ) , respectively.  We  shall  see  that  wave  propagation 
in  the  basic  atmosphere  may  be  completely  characterized  by  the  vertical  variation  of  the  sound  speed  c(z),  the 
Alfven  velocity  i a(j),  and  the  local  density  scale  height  H(z),  defined  by 


The  unperturbed  atmosphere  is  taken  to  be  in  static  equilibrium: 


If  the  magnetic  fi;ld  is  a function  of  height  then  it  has  a role  in  the  basic  equilibrium  of  the  atmosphere.  We 
consider  only  stable  atmospheres,  which  requires  that 

‘JTo  _ (dM  ]_±  (Bl\ 

dz  \ dz  /s  p0R  dz  \ %ir  } 

where  ( dTu;dz)s  is  the  adiabatic  temperature  gradient  and  R is  the  gas  constant  (see  Thomas  and  Nyc  1975  for  a 
recent  discussion) 

Consider  small  adiabatic  perturbations  of  the  equilibrium  atmosphere,  letting  p,  p,  «,  and  B denote  the  per- 
turbations in  pressure,  density,  velocity,  and  magnetic  field,  respectively.  Then,  the  basic  linearized  equations  of 


continuity,  momentum,  energy,  and  induction  are 

| + V.U«)  = o,  (3) 

Po  Tt  + ~ pg  ~ ^ * Bo)  * B + (V  X B)  X B0]  = 0 , (4) 

^ + u.VPo  = c*(j?  + u VPoy  ' ( (5) 

< B 

— - V x (u  x B0)  = 0 . (6) 
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After  taking  the  time  derivative  of  the  momentum  equation  (4),  we  may  eliminate  the  perturbation  quantities 
p,  p,  and  B by  using  equations  (3),  (5).  and  (6)  This  leaves  a single  vector  equation  for  the  velocity  perturbation 
u = (tt,  t\  H'l: 

p0Wi  + V[Po’‘'(^  + g~~7r)  “ + [(w  V)Po  + p0V  u\g 

- x B°)  x[V  x («  x if0)]  - B0  x (V  x [V  x (u  x £„)]}]]  = 0.  (7) 

Next,  we  assume  that  the  perturbation  velocity  has  the  form  u = u<ixpi(k  r - u>t),  with  u = u(z)  = [ii(z),v(z),  w(:)] 
and  k r = kxx  + kyy.  Then,  using  the  definition  of  eA2,  the  three  components  of  the  momentum  equation  (7) 
become 

(w2  - c2k2)u  - c2kxkvv  - — c2  >v  = 0 , (8) 

—c2kxkyii  + [ui2  - c2k2  - vA2(k2  + k2)]v  - ik^g  - (e2  + eA2)^  w = 0 , (9) 


Idc2  c 2 2 d\  • , •/  f d 

'*x\dz  H + S + C dz)U  + ,k\dz 


(c2  + i-A2)  - 


2.  _ (c2  + i\2) 


+ g + (c2  + eA2)  jz  v 


+ (c2  + t’A2) 


d2w  d 
+ [dz 


(e2  + vA2)  - 


(c2  + i>A2)] 


+ (w2  - vA2kx2)w  = 0 . (10) 


These  are  the  linearized  perturbation  equations.  They  give  important  information  about  particle  motions  for  various 
modes  of  propagation  in  the  atmosphere. 

The  horizontal  components  of  the  perturbation  velocity  can  be  eliminated  from  the  system  of  equations  (8)— ( 1 0) 
to  yield  a single  equation  for  the  vertical  velocity  w.  The  resulting  equation  is 

S + A(z)^  + B(z)w  = 0,  (11) 


where  the  coefficient  A (z)  and  B(z)  are  given  by 

A{Z)  = “77  + ~ ^2^2)2  £ + i[-^  + (o>2  - c2kx2)(w2  - vA2k2) 


B(z)  = ^ j u>6  - [(c2  + rA2)(kx2  + ky2)  + vA2k2)^ 


+ v2kx\k2  + k2){2c2  + rA2)  - g(k2  + *,•)(*  *AV  a,2 

- V\k2(k2  + kv2)cWk2  -*(*-■£)  ~§  ("2  ~ vA2kx2)2w2(kx2  + kS)d-£  - ~ kv2g • 


Here,  D and  E are  given  by 


D = (<u2  — cA2^2)[a>2(c2  + t>A2)  - c2l’A2/cx2] 


£ = cu-  - (*x2  + &v2)[u>2(c2  + vA2)  - c2rA3fcx2] . (15) 

Equation  (II)  is  a general  equation  for  the  vertical  component  of  velocity  for  a perturbation  propagating  in 
an  arbitrary  direction,  in  an  atmosphere  with  a horizontal  magnetic  field  and  arbitrary  vertical  distributions  of  c2, 
vA  , and  H.  This  propagation  equation  has  been  given  previously  (Nye  and  Thomas  1974u)  in  the  case  ky  = 0. 
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Now  consider  the  case  where  the  undisturbed  temperature  and  magnetic  field  are  constant  with  height.  Since  the 
magnetic  field  is  uniform,  it  has  no  effect  on  the  hydrostatic  equilibrium  of  the  atmosphere,  and  the  equilibrium 
pressure  and  density  both  decrease  exponentially  with  height.  The  sound  speed  and  the  density  scale  height  are 
both  constant,  with  values  determined  by  the  temperature  of  the  atmosphere.  The  Alfven  velocity  increases 
exponentially  with  height  due  to  the  decreasing  density.  The  sound  speed,  density  scale  height,  density,  and  Alfven 
velocity  are  given  by 

c = ( yRT0 )1/2  = const. , (16) 


and 


H = — = const. , 
g 

Po(:)  = Poo?  *,H  . 
Da(z)  = I'oe2'2" , 


(17) 

(18) 
(19) 

where  p00  and  i>0  are  the  values  of  the  undisturbed  density  and  Alfven  velocity  at  z = 0. 

The  nondimensional  parameter/?2  = v02/c2  is  introduced  as  a measure  of  the  relative  importance  of  the  restoring 
forces  due  to  the  magnetic  field  and  to  compressibility  at  the  point  z = 0.  For  values  of  /32  < 1 , there  is  a region  of 
the  atmosphere  above  z = 0 where  compressibility  has  more  importance  as  a restoring  force  than  does  the  magnetic 
force,  but  in  any  case  the  magnetic  field  always  becomes  dominant  as  z becomes  large. 

In  the  remainder  of  this  paper  we  shall  consider  only  waves  whose  horizontal  component  of  propagation  is 
parallel  to  the  magnetic  field  (k„  = 0).  Using  (I6)-(19),  we  may  write  (1 1)  as 


[c2w2  + (co2  — c2kx2)  v02e*IH] 


d2w  c2tu2  dw 


dz2 


H dz  + 


(w2  - c2k2)( cu2  - i-o 2e*'Hk2)  - g^g  - Jjjkx2  w = 0. 


We  define  the  nondimensional  frequency  £2  and  the  nondimensional  horizontal  wavenumber  A'  by 

Q = Hoj/c  , K = Hkx  . 

By  transforming  the  dependent  and  independent  variables  according  to 

Q.2 

Lf/  — / H y — fj—Z/H 

’ /?2(A2  - ii2) 


we  may  put  equation  (20)  in  the  dimensionless  form 

d2W 


,dW 


with 


41  ~ X)l&  + [C-  ^ + B + 1)*]-^  - ABW  = U, 


A + B = C = 2A  + 1 , AB  = a2  + K + § 


(20) 


(21) 


(22) 


(23) 


(24) 


Equation  (23)  is  the  standard  form  of  the  hypergeometric  differential  equation.1  The  solutions  of  this  equation 
may  be  expressed  in  terms  of  hypergeometric  functions,  given  for  |x|  < 1 by 


F(A,  B;  C;  x)  = 


r(C)  v r(/f  + n)V(B  + n)  x* 
r(A)l\B)  „40  T(C  + n)  nl 


(25) 


The  general  solution  of  equation  (23)  may  be  written,  for  |x|  < 1,  in  terms  of  the  original  variables  z and  w(z) 


as 


w(z)  = D1e-*KIHF 


A,  B\  C; 


a2 


+ D2e*KIH 


P\K2  - U2) 

/32(A2  - O2) 


£22 


' 2K 

F A - 


C + 1,5-  C+  1 ; 2 - C; 


n2 


P2(k2  - a2) 


(26) 


where  D,  and  D2  are  arbitrary  constants. 


1 The  hypergeometr x nature  of  the  wave  equation  in  this  case  was  noted  by  us  earlier  (Nye  and  Thomas  19746),  and  also  in- 
dependently by  Adam  < 1975). 
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IV.  E1GENMODES  FOR  A RIGID  LOWER  BOUNDARY 

We  now  examine  modes  of  propagation  in  the  relatively  simple  case  of  an  isothermal  atmosphere  with  a uniform 
horizontal  magnetic  field,  bounded  from  below  by  a rigid  wall  at  z = 0.  The  general  solution  (26)  is  subject  to 
boundary  conditions  at  z = 0 [x  = ft2/j32(A'2  - ft2)]  and  at  z = oo  (x  = 0). 

As  a condition  at  z = oo,  we  require  that  the  total  energy  of  the  perturbation  be  integrable  over  0 < z < oo. 
The  magnetic  energy  of  the  perturbation  is  proportional  to  the  square  of  the  velocity.  Since 

lim  F(a,  P;  y;  x)  - 1 , (27) 

x-0  7 

we  see  from  equation  (26)  that  we  must  take  D2  = 0 . 

The  second  boundary  condition  is  that  the  vertical  velocity  vanish  at  the  rigid  wall,  i.e.,  w = 0atz  = 0ie  at 
x = Q2/P2(K2  - ft2).  Provided 


ft2 

P2(K2  - ft2)  I 41  1 ’ (28) 

we  may  apply  this  condition  directly  to  equation  (26)  with  D2  = 0 to  obtain  the  dispersion  relation 


fl;C; 


= 0. 


(29) 


If,  however,  (28)  is  not  satisfied,  then  other  representations  of  the  general  solution  (26),  valid  for  \x\  > 1,  must  be 
used  in  order  to  apply  the  boundary  condition  at  z = 0.  For 


1 

P \K2  - ft2)  ~ ’ (3°) 

the  interval  0 < x < Ll2/p2(K2  - ft2)  contains  the  regular  singular  point  at  x = 1,  and  no  solution  that  satisfies 
the  boundary  conditions  and  is  also  regular  at  x = 1 is  to  be  expected.  For  the  range 

ft2  . , 

P2(K2  - ft2)  " 1 ’ <31) 


we  may  use  analytic  continuation  to  extend  the  general  solution  (26).  The  analytic  continuation  of  (26)  with  D2  = 0, 
valid  for  x < — 1,  is  given  by 


w(z)  = ZV'**'"- 


r fr(or(5  — a) 

\ [iw 


p\n2  - k2) 


ft2 


e*AmF\A,  1 — B\  1 — B + A; 


_ p2(K2  - ft2) 


ft2 


gZIH 


r(C)F(A  - B) 

[r(A)}2 


P2(Q2  - K2) 
ft2 


^b,hf 


Thus,  for  the  rtnge  of  parameters  (31),  the  dispersion  relation  is  given  by 


(32) 


T(g  - A) 

[F(g)]2 


^2(ft2  - a:2)] 


ft2 


V[/(,  1 - B;  1 - B + A ; Q8)| 


r (a  - b) 

[r(A)]2 


jS2(ft2  - A2)] 
ft2 


= 0 . (33) 


The  dispersioi  relation  (eq.  [29]  or  [33])  has  been  evaluated  for  various  values  of  the  nondimensional  parameter 
p2.  Examples  are  plotted  in  Figures  1 and  2 (also  sec  Fig.  4).  The  curves  in  these  figures  represent  well-defined 
eigenmodes  of  wave  propagation  in  the  atmosphere.  These  curves  represent  trapped  waves  propagating  horizont- 
ally in  the  wave  guide  formed  by  the  solid  boundary  below  and  the  exponentially  increasing  Alfven  velocity  above 
There  are  several  ways  of  interpreting  the  effect  of  changing  P2  on  the  dispers'ion  relation.  First,  different  values 
of  P2  can  be  taken  to  represent  the  same  magnetic  field  strength  and  the  same  density  at  z = 0,  but  different  atmos- 
pheric temperatures.  A second  interpretation  is  that  different  values  of p2  represent  the  same  magnetic  field  strength 
and  the  same  ten  peraturc,  but  different  densities.  This  is  equivalent  to  placing  the  solid  lower  boundary  at  succes- 
sively higher  leve*s  in  the  atmosphere  corresponding  to  larger  values  of  /?2.  At  each  higher  level  the  magnetic  restor- 
ing force  becomes  more  important  due  to  the  decreased  density,  while  the  compressible  restoring  force  remains  the 
same. 
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Fig.  1. — Dispersion  diagrams  (nondimensional  frequency  versus  nondimcnsional  horizontal  wavenumber)  for  an  isothermal 
atmosphere  with  a uniform  horizontal  magnetic  field  and  a solid  lower  boundary,  with  y = 5/3,  /S2  = I0's,  and  fp  = 10'4.  The 
curves  represent  eigenmodes,  and  the  crosses  indicate  the  computed  points.  This  figure  should  not  be  confused  with  a diagnostic 
diagram  for  an  atmosphere  with  constant  parameters  (e.g.,  Yu  1965). 

Fig.  2. — Same  as  Fig.  1,  but  with  /32  = 0.5. 


The  third  interpretation  is  to  consider  changes  in  )32  to  be  due  to  changes  in  the  magnetic  field  strength,  with 
fixed  values  of  temperature  and  density.  Since  the  sound  speed  and  the  density  scale  height  then  do  not  change,  the 
scales  for  the  frequency  and  horizontal  wavenumber  are  the  same  in  each  case  and  the  dispersion  diagrams  can  be 
compared  directly.  From  Figures  1 and  2 it  can  be  seen  that  increasing  the  magnetic  field  strength  (increasing  /l2) 
increases  the  cutoff  frequency.  As  ft2  increases,  the  slope  of  the  dispersion  curves,  and  hence  the  group  velocity, 
also  increases. 

Lowering  the  value  of  y to  represent  crudely  the  effect  of  radiative  transfer  has  little  effect  on  the  nondimensional 
dispersion  diagrams.  However,  the  frequency  scaling  depends  on  y as  o>  ~ (y)l'2D,  while  the  wavenumber  scaling 
is  independent  of  y.  Therefore,  for  lower  y (lower  sound  speed)  the  wave  oscillates  less  rapidly  and  the  phase  and 
group  velocities  are  correspondingly  lower. 

The  vertical  velocity  of  the  disturbance  (eq.  [26]  or  [32])  can  be  calculated  as  a function  of  height  for  any  point 
on  a dispersion  curve.  Figure  3 compares  the  lowest  mode  of  oscillation  for  the  same  horizontal  wavenumber  but 
different  values  of/32  (i.e.,  different  magnetic  field  strengths),  and  shows  that  for  increasing  magnetic  field  strength, 
the  wave  oscillates  more  rapidly  and  is  trapped  at  lower  levels  in  the  atmosphere.  We  now  discuss  the  eigenmodes 
given  by  the  dispersion  relation  (29)  or  (33)  in  relation  to  an  observed  solar  oscillation. 

V.  APPLICATION  TO  FLARE-INDUCED  CORONAL  WAVES 

On  1963  September  20,  Moreton  and  Ramsey  (Moreton  1965)  observed  a chromospheric  disturbance,  apparently 
caused  by  the  flasfi  phase  of  a flare,  propagate  at  a nearly  constant  velocity  of  750  kms"1  for  several  hundred 
thousand  kilometers  across  the  solar  disk.  Many  other  flare-induced  disturbances  have  been  reported  (Moreton 
1960;  Athay  and  Moreton  1961 ; Dodson  and  Hedeman  1968),  and  the  propagation  velocity  is  usually  on  the  order 
of  1000  km  s'1.  Dodson  and  Hedeman  (1968)  report  that  the  width  of  the  disturbance  created  by  the  proton  flare 
of  1966  August  28  was  greater  than  100,000  km. 

These  disturbances  could  not  have  been  propagating  solely  in  the  chromosphere,  since  in  the  chromosphere  the 
sound  speed  is  only  of  the  order  of  20  km  s'1  and  the  Alfven  velocity  is  only  of  the  order  of  50  km  s-1.  Thus  a 
purely  chromospheric  disturbance  would  have  created  a shock  wave  and  been  rapidly  dissipated.  In  the  corona, 
however,  both  the  Alfv£n  velocity  and  the  sound  speed  are  an  order  of  magnitude  higher  than  in  the  chromosphere 
due  to  the  increased  temperature  and  decreased  density.  It  has  been  proposed  (Meyer  1968;  Uchida  1968,  1970, 
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Fig.  3.— Direct  comparison  of  the  first  mode  of  oscillation  of  the  atmosphere  for  the  same  nondimensiona!  horizontal  wave- 
number  (A  = 0.25)  but  for  two  values  of  fP  (fi2  = 10'3,  j33  = 10'1). 

1974;  Uchida  et  al.  1973)  that  the  disturbance  is  a magnetohydrodynamic  wave  propagating  in  the  low  corona  and 
that  the  motion  of  this  wave  at  the  corona-chromosphere  transition  region  is  what  is  actually  observed.  There  is, 
however,  no  general  agreement  as  to  the  wavelength  of  the  disturbance  or  even  whether  the  observed  disturbance 
is  a single  wave  or  a wave  packet. 

Meyer  ( 1 968)  studied  the  propagation  of  the  magnetoacoustic  fast  mode  in  an  isothermal  corona  permeated 
by  a uniform  vertical  magnetic  field,  with  a rigid  lower  boundary  representing  the  chromosphere-corona  transition 
region.  He  found  eigenmodes  with  nearly  constant  horizontal  group  velocity.  Equating  the  group  velocity  to  the 
observed  propagation  velocity,  Meyer  found  that  for  a horizontal  wavelength  on  the  order  of  100,000  km,  the 
coronal  magnetic  field  must  be  approximately  6 gauss,  a reasonable  average  value. 

Uchida,  in  a series  of  papers  (Uchida  1968,  1970,  1974;  Uchida  et  al.  1973),  studied  the  propagation  of  short- 
wavelength  ( ~ 5000  km)  disturbances  in  various  realistic  coronal  models.  Using  a ray-tracing  technique,  he  obtained 
horizontal  and  vertical  refractions  in  close  agreement  with  the  observed  waves. 

Although  the  magnetic  field  structure  of  the  corona  is  quite  complicated,  the  field  changes  fairly  slowly  and  there 
are  probably  regions  of  nearly  uniform  field  with  almost  any  orientation.  We  study  the  case  of  a uniform  horizontal 
field  in  connection  with  the  coronal  wave  problem  only  as  a means  of  understanding  the  mechanism  of  wave 
propagation  for  waves  of  arbitrary  wavelength.  Our  model  supplements  Meyer’s  (1968)  work  by  considering  the 
case  of  a uniform  horizontal  magnetic  field,  and  by  including  effects  of  gravity  and  stratification.  As  in  Meyer's 
model,  we  use  the  rigid  lower  boundary  to  represent  upward  reflection  from  the  chromosphere-corona  transition 
layer. 

We  have  evaluated  our  solution  for  a temperature  of  1.6  x 10eKand/32=  10.whichisfairlytypicaloftheba.se 
of  the  corona.  These  values  correspond  to  a sound  speed  of  180  km  s' 1 and  a density  scale  height  of  approximately 
71,000  km.  Figure  4 shows  the  dimensional  dispersion  relation  for  these  parameters.  We  see  that  the  dispersion 
curves  are  nearly  straight,  which  means  that  these  modes  have  very  little  dispersion  and  will  propagate  for  great 
distances  with  little  change  in  character.  The  first  three  modes  have  been  calculated  for  a wavelength  of  75,000  km 
in  Figure  5.  The  first  mode  has  nearly  zero  vertical  velocity  above  two  scale  heights  and  is  therefore  trapped  in 
the  low  corona. 

The  phase  velocity  and  group  velocity  of  the  first  mode  have  been  plotted  as  a function  of  horizontal  wavenumber 
in  Figure  6.  For  any  wavelength  of  100,000  km  or  less,  the  group  velocity  is  nearly  constant  at  about  610  km  s'1. 
Since  the  energy  of  a disturbance  propagates  at  the  group  velocity,  it  is  not  important  which  specific  wavelength, 
or  spectrum  of  wavelengths,  receives  energy  from  the  flare.  The  energy  at  all  wavelengths  will  propagate  together 
as  a wave  packet  near  the  lower  coronal  boundary. 
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Fig.  4. -Dimensional  dispersion  diagram  for  fP  = 10,  y = 5/3,  and  T0  = 1.6  x 106  K.  The  curves  represent  eigenmodes  of  the 
corona  which  are  trapped  by  the  increasing  Alfven  velocity  with  height. 


The  present  model  is  not  proposed  as  a realistic  model  of  the  solar  corona,  although  it  may  be  fairly  accurate 
over  certain  regions.  No  attempt  has  been  made  to  include  the  elfects  of  horizontal  variations.  The  value  of  the 
model  is  that  a mechanism  lor  wave  propagation  can  be  studied  for  arbitrary  wavelengths.  These  results  close  the 
gap  between  the  short-wavelength  ray-tracing  theory  and  the  long  wavelength,  vertical  field  case.  We  show  that 
the  question  of  wavelength  is  not  particularly  important  since  the  group  velocity  of  the  trapped  modes  is  essentially 
independent  of  wavelength.  For  the  relatively  large  value  of  #2  ( = 10),  the  wave  modes  are  basically  the  magneto- 
acoustic  last  maces  (studied  by  Meyer  and  by  Uchida)  modified  by  gravity.  For  an  inclined  magnetic  field,  there 
must  also  be  trapped  modes  of  propagation  involving  a coupling  of  the  present  modes  and  tne  type  of  mode 
studied  by  Meyer  for  a vertical  field.  The  present  results,  taken  with  "those  of  Meyer  and  the  work  of  Uchida,  present 
a consistent  pictu-e  of  flare-induced  coronal  waves  as  guided  magneto-atmospheric  waves. 


Fig.  5.  First  three  modes  of  oscillation  of  the  model  corona  for  p2  = 10,  y = 5/3.  The  vertical  velocities  have  been  normalized 
to  maximum  value  unity. 
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Fig.  6. — Phase  velocity  ep  and  group  velocity  v,  of  the  first  mode  of  coronal  oscillation  plotted  as  a function  of  horizontal 
wavenumber. 
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ABSTRACT 

A simple  two-layer  model  of  a sunspot  penumbr*  is  used  to  study  the  mode  of  running  pen- 
umbraf  waves.  Exact  solutions  of  the  linearized  wave  equation,  not  limited  to  the  small-wavelength 
approximation,  are  employed  in  each  layer.  The  lowest  “plus”  eigenmode  of  magneto-atmospheric 
waves  in  the  model  penumbra  is  in  good  agreement  with  observations  of  running  penumbral  waves. 
The  results  indicate  that  running  penumbral  waves  should  be  observable  in  a photospheric  spectral 
line. 

Subject  headings:  hydromagnetics  — Sun:  atmospheric  motions  — Sun:  magnetic  fields  — 

Sun:  sunspots 


I.  INTRODUCTION 

In  Paper  1 of  this  series  (Nye  and  Thomas  1976)  we  presented  an  exact  analytical  solution  for  magneto-atmos- 
pheric waves  in  the  case  of  an  isothermal  atmosphere  with  a uniform  horizontal  magnetic  field.  In  the  present  paper 
we  apply  this  solution  to  a simple  two-layer  model  penumbra  in  order  to  study  the  mode  of  running  penumbral 
waves. 

Running  penumbral  waves  (Zirin  and  Stein  1972;  Giovanelli  1972,  1974;  Moore  and  Tang  1975)  are  good 
examples  of  magneto-atmospheric  waves.  These  waves  propagate  radially  outward  across  sunspot  penumbrae, 
with  predominantly  vertical  motions  in  Ha.  The  observed  range  of  frequency  and  propagation  speed  is  fairly  well 
established  (see  discussion  in  § IV). 

Moore  (1973)  has  concluded  that  the  source  of  excitation  of  the  penumbral  waves  is  overstable  convection  in  the 
low  umbra.  In  an  earlier  paper  (Nye  and  Thomas  1974  [NT])  we  studied  the  mode  of  propagation  of  penumbral 
waves  on  the  basis  of  a piecewise  linear  model  of  the  vertical  structure  of  a typical  sunspot  penumbra.  We  found 
the  penumbr.il  waves  to  be  magneto-atmospheric  waves  (of  the  “plus”  type)  that  are  vertically  trapped  at  photo- 
spheric  levels.  This  trapping  is  primarily  due  to  the  increasing  sound  speed  with  depth  into  the  convection  zone 
and  the  increasing  Alfven  velocity  with  height  into  the  chromosphere. 

Here  we  ex  end  our  earlier  work  by  computing  actual  eigenmodes  of  propagation  for  a somewhat  simpler  model 
penumbra,  which  nevertheless  retains  the  essential  features.  The  properties  of  the  lowest  mode  of  propagation  of 
the  model  penumbra  turn  out  to  be  in  good  agreement  with  observations,  and  give  some  useful  clues  for  further 
observation  of  running  penumbral  waves. 

II.  THE  TWO-LAYER  PENUMBRAL  MODEL 

The  entire  penumbral  mode!  consists  of  a compressible,  inviscid,  perfectly  conducting,  stratified  perfect  gas 
subject  to  a constant  acceleration  of  gravity  g ( = 0.274  km  s'2)  in  the  negative  c-dircction.  The  upper  layer  is 
isothermal  and  is  permeated  by  a uniform  horizontal  magnetic  field,  which  yields  an  Alfven  velocity  that  increases 
exponentially  with  height  due  to  the  decreasing  density.  An  exact  solution  of  the  linearized  propagation  equation 
for  this  case  was  given  in  Paper  I.  i 

This  upper  layer  is  a suitable  model  of  the  penumbrai  photosphere  and  chromosphere,  where  observed  penumbral 
magnetic  field;  are  very  nearly  horizontal  (Nishi  and  Makita  1973)  and  decrease  slowly  with  height  (Bray  and 
Loughhead  1954).  The  scale  height  for  variation  of  the  magnetic  field  is  very  large  compared  to  the  density  scale 
height,  so  the  assumption  of  the  uniform  horizontal  magnetic  fielc  is  reasonable.  Our  earlier  calculations  (NT) 
showed  that  running  penumbral  waves  are  trapped  at  photospheric  levels,  so  that  the  increasing  sound  speed  in 
the  upper  chromosphere  has  little  effect  on  the  trapping.  Taking  the  upper  layer  to  be  isothermal  is  therefore  also  a 
reasonable  assumption. 

The  vertical  distributions  of  the  sound  speed  and  Alfven  velocity  for  ti  e two-layer  model  penumbra  are  shown 
in  Figure  1.  Subscripts  1 and  2 denote  quantities  in  the  upper  and  lower  layers,  respectively,  and  the  subscript  0 


* Also  C.  E.  Kenneth  Mces  Observatory. 
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P‘str'bution  of  c3,a"d  !i*2  *ith  height  z in  the  two-layer  penumbral  model.  The  upper  layer  (1)  is  isothermal  with 
uniform  horizontal  magnetic  field.  The  lower  layer  (2)  has  an  adiabatic  temperature  gradient  and  no  magnetic  field. 


refers  to  quantities  evaluated  at  z 
layer  are  given  by 


- 0.  The  sound  speed,  density  scale  height,  and  Alfven  velocity  in  the  upper 


and 


Cia  = yiRTl  = const. , 

Hi  = — = const. , 
yig 


(1) 

(2) 


va2(z)  = v02  exp  (z/Hi) . 


V) 


The  lower  layer  of  the  penumbral  model  (layer  2,  figure  1)  is  adiabatic  with  no  magnetic  field.  The  temperature 
decreases  with  height  (increases  with  depth)  at  the  adiabatic  lapse  rate,  (dT/dz)s  = -g/cp,  and  thus  this  layer  is 
neutrally  stable.  The  actual  temperature  distribution  in  the  convection  zone  below  a penumbra  is  probably  very 
nearly  adiabatic,  except  for  a thin  superadiabatic  layer  just  beneath  the  photosphere  that  we  neglect  here  There  is 
no  magnetic  field  in  this  layer  since  we  assume  that  the  penumbral  magnetic  field  lies  over  the  convection  zone. 
I he  sound  speed  squared  and  local  density  scale  height  each  increases  linearly  with  depth  in  the  lower  layer  their 
functional  forms  being 


and 


c*2(z)  = c202  - g(y2  - l)z 


(4) 


H2(z)  — H2 o ~ (y2  — 1)2  • 

The  corresponding  density  distribution  is 


(5) 


p2(z)  = f>2o[l  - (y2  - l)z///20]1,<’,»-1) . (6) 

At  the  interface  between  the  two  layers  (2  = 0),  we  require  the  undisturbed  density  to  be  continuous  to  avoid 
introducing  interfacial  gravity  waves  and  wave  reflections;  therefore,  Pl0  = p20.  In  the  unperturbed  penumbra, 
there  must  be  pressure  equilibrium  at  the  interface;  that  is,  the  gas  pressure  in  the  lower  layer  at  2 = 0 must  equal 
the  sum  of  the  gas  pressure  and  the  magnetic  pressure  in  the  upper  layer  at  z = 0.  Therefore,  the  gas  pressure  is 
greater  in  layer  2 than  in  layer  1 ; and  since  the  density  is  continuous  across  the  interface,  the  temperature  is  greater 
in  layer  2.  This  may  be  expressed  in  terms  of  the  sound  speeds  and  the  Alfven  velocity  at  z = 0 as 


r 2 - l2  r 2 4.  ' 2 2 

c20  - ~ ci  + y i>o  - 


(7) 


We  now  turn  to  the  problem  of  computing  eigenmodes  of  magneto-atmospheric  waves  in  this  penumbral  model. 


111.  ANALYSIS 


Consider  first  the  behavior  of  small  adiabatic  perturbations  in  the  lower  layer  (layer  2),  an  adiabatic  atmosphere 
without  magnetic  held.  Leibachcr  ( 1971)  solved  this  problem  in  his  study  of  oscillations  of  the  quiet  photosphere. 
Here  we  take  a slightly  different  approach  from  his,  using  other  transformations  which  yield  a different  form  of  the 
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propagation  equation.  For  vanishing  magnetic  field  (B0  = B 0),  the  vector  equation  for  the  perturbation 
velocity  (eq.  [7],  Paper  1)  becomes 


<-u2 

<Jt2 


= + (72  - 1 Hg  + ^(«2  g)  • 


(8) 


where  <t>  = V u2. 

We  assume  that  the  perturbation  velocity  has  tije  form  u2  = u^(;)e.\p  [i(kxx  — to/)],  with  propagation  in  the 
^-direction  (A„  0).  This  implies  that  f2  = i)  and  u„.  -=  i/Jr)  - [u2(z),  0,  h’2(j)].  All  other  perturbation  quantities 

are  represented  in  a similar  manner,  with  a caret  denoting  the  r-dependent  amplitude  in  each  case.  From  the  two 
components  of  equation  (8)  and  the  definition  of  we  obtain  the  following  relation: 


(z)  = 


Kb*  ~ c22kx2lw2]<j,  - c22d<f>/dz 


(9) 


<o2  - g2kx2/< O2 

The  pressure  perturbation  can  also  be  written  in  terms  of  <f>,  using  the  continuity  and  energy  relations  (see  Paper  I): 


'P2<2 


(w2  - g2y2/c22)4>  + gd<f>/dz 

w2  - g 2kx2/oj2 


(10) 


Upon  substitution  of  equation  (9)  into  the  e-component  of  equation  (8),  we  obtain  the  following  second-order 
differential  equation  for  if>: 


_ .wA  ‘Hi 

dz2  \cf  dz  c22 ) dz 


i-  2 

2 

(2 


+ s4h<r,-  0 + ^^' 


CL>  C2 


<j>  = 0. 


(11) 


The  nondimensional  frequency,  horizontal  wavenumber,  and  depth,  based  on  the  values  of  sound  speed  and 
density  scale  height  in  layer  2 at  z = 0,  are  defined  by 


12-2  = 7/2001/020  , *2  3 // 2i>k x , and  z = z/ H 2( 

By  transforming  the  independent  and  dependent  variables  according  to 

2*2 


Y = 


V -2  ~ I 


~ 2*2-% 


= ey'2<P , 


equation  (11)  assumes  the  form 


where 


, d2^ 

dY2 


d<l> 


+ (b-  = 0, 


(2 y2  - 1) 

2(g2  - 1)  2*2(72-0 


b = 


(2y,  ~ 0 
(72  - 0 


(12) 


(13) 


(14) 


(15) 


Equation  (14)  is  the  standard  form  of  Kummer’s  equation  (see  Abramowitz  and  Ste gun  1964).  The  solutions  of 
this  equation  are  given  in  terms  of  Kummer’s  functions, 


r(fl)B4o  1 \b  + n)  n\ 

The  general  solution  of  equation  (14)  is 

•K  Y)  = b\  Y)  + D<U(a . b;  Y) , 


(16) 


(17) 


where  D-,  and  D,  are  arbitrary  constants  and  (/(a,  b\  Y)  can  be  written  in  terms  of  Rummer’s  functions.  To  insure 
finite  total  perturbation  energy,  we  must  require  that  the  vertical  velocity  of  the  perturbation  vanish  as  z -*■  —a 
or  as  Y * +00.  This  in  turn  requires  that  /):l  = 0.  The  solution  for  <f>  as  a function  of  z in  the  lower  layer  is  then 


m = £>, 


sin  -rrb 


exp  i - 


*2 


(72  - 1) 


— KoZ 


M[a,  6;  2KJ(y2  - 1)  - 2 *25] 
l’(l  -ha-  b)V(b) 


2*2 


1 (72  I) 


- 2*27 


■-'>  M[\  + a - b.  2 - 6;  2 *2/(72  - 1)  - 2K2z]'\ 

I’(fl)r(2  - b)  f 


(18) 


33 


No.  2,  1976 


SOLAR  MAGNETO-ATMOSPHERIC  WAVES 


585 


The  vertical  velocity  and  pressure  perturbation  are  given  in  terms  of  by  equations  (9)  and  (10). 

The  term  of  the  solution  in  the  upper  layer  has  been  given  in  Paper  1,  and  we  shall  not  repeat  the  analysis  here. 
The  general  form  tor  the  vertical  velocity  in  the  upper  layer  that  gives  finite  total  perturbation  energy  is  given  in 
terms  of  hypergeometric  functions  by  either 

r Qi 

I A,  B,  C , a2l  2 ( \ 2~\  ®XP  ( — z///j) 


for 

or  else  by 
w,(z)  = D 


*i{z)  = exp  (-zKJHl)F 

«i2 


0W  - <V) 

fi\KS  - Q,2)  1 < eXp  ’ 


09) 

(20) 


fr(C)T(B  - A) 

m2  - KS)} 

1 [IX*)]2 

L ‘V 

x F 

+ 


|/f,  1 - B;1  - 

T(C)r (A  - B) 

[IW 


fi  + ,WMV)cxp(2W 

nns  - 


X F 


2 _ K 

7^2 — —J  exp  (zB/Hi) 
B,l  - A ;1  - A + B ; 


.^W-JV) 

o? 


exp  (z///t) 


for 


Here 


ns 


> exp  (z/HS) . 


I W - <V) 

A + B = C = 2K,  4-  1 , 
AB  . 0,*  + AT,  + (^)  | 


(21) 

(22) 

(23) 

(24) 

(25) 


jS2  = VoVcS  , 

and  Q!  and  are  nondimensional  frequency  and  wavenumber  defined  as  in  (12),  except  sealed  with  0 and  // 
The  pressure  perturbation  in  the  upper  layer  consists  of  the  sum  of  the  gas  pressure  perturbation  and  the  mag- 
netic pressure  perturbation.  The  gas  pressure  perturbation  can  be  expressed  in  terms  of  the  vertical  velocity  and 
its  derivative  as 


Pi  = -ipicSu) 


Idwy  g „ \ 


V - cSkS) . 


(26) 


The  magnetic  pressure  perturbation  pm  is  found  from  the  linearization  of 

_ , , (tf0  + B)-(B0  + B) 

Pm  ' Pm  » 


where  p„  is  the  unperturbed  magnetic  pressure.  The  components  of  the  perturbed  magnetic  field  are  determined 
by  the  linearized  induction  equation  (eq.  [4]  of  Paper  I).  We  find 


Pm  — 


. va2Pi  dw1 

w dz 


(27) 


The  total  perturbed  pressure  pTl  in  layer  1 is  expressed  in  terms  of  the  density  and  the  Alfven  velocity  as 

Ptx  = ~ -j  l(cS  + vK2W  - cSvSk*)  — - w2gw^j(w2  - cSkx 2) . (28) 

\ 

IV.  EIGENMODES  AND  RUNNING  PENUMBRAL  WAVES 

We  now  have  expressions  for  the  vertical  velocity  and  the  pressure  perturbation  in  each  layer  of  the  penumbral 
model,  such  that  the  total  perturbation  energy  is  finite.  The  remaining  conditions  are  the  matching  of  the  vertical 
velocity  and  the  perturbed  pressure  at  the  interface  z = 0.  1 
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The  scaling  of  frequency  and  wavenumber  was  done  separately  for  each  layer  in 
pagation  equation  as  much  as  possible  in  each  case.  In  matching  across  the  interface, 
lations  between  parameters  in  the  two  layers: 

order  to  simplify  the  pro- 
we  need  the  following  re- 

Hia  = y2(i  + i^yJZ-A  , 

(29) 

U22  = 7W1  + WnW  - 

(30) 

= y22(  1 + WnW . 

and  equation  (7). 

Continuity  of  the  vertical  velocity  across  : = 0 requires,  after  normalization,  that 

(31) 

HT(0)  = tv2(0)  = 1 . 

(32) 

This  condition  fixes  the  values  of  the  coefficients  Z),  in  equation  (19)  or  (21)  and  Z)4  in  equation  (18).  The  remaining 
condition,  the  continuity  of  the  perturbed  pressure,  requires  that  we  equate  (10)  awd  (28).  This  then  leads  to  the 
nondimensional  condition 


[(i  + pw  - PKMd&jdzju  - \_ 
iii2  - KS  ) Yi 

where  Equation  (33)  is  only  satisfied  by  particular  values  of  frequency  and  wavenumber,  and  gives  the 

dispersion  relation  for  eigenmodes  of  oscillation  in  the  penumbral  model. 

In  order  to  evaluate  the  dispersion  relation  (33).  the  free  parameters  /32,  y,  and  y2  must  be  specified,  which  then 
effectively  determines  the  properties  of  the  model.  Although  it  is  possible  to  evaluate  (33)  for  different  values  of 
y in  each  layer  (for  example,  a lower  value  of  yl  could  be  taken  to  represent  radiative  transfer  in  the  upper  layer), 
we  chose  the  usual  value  of  5/3  for  both  layers.  The  dispersion  relation  (33)  was  solved  numerically  by  inserting 
values  of  A'i  and  then  computing  and  comparing  the  two  sides  of  (33)  for  small  increments  in  . 

Figure  2 shows  the  first  several  eigenmodes  of  the  two-layer  penumbral  model  for/22  = 0.5  and  y;  = y2  = 5/3. 
The  value  of  /32  was  chosen  to  represent  a typical  penumbra  and  is  slightly  less  than  the  value  of  /32  at  : = 6 in  our 
earlier  penumbral  model  (NT).  Here  we  have  classified  the  eigenmodes  as  "plus"  or  “minus”  modes,  following 
the  terminology  used  in  the  case  of  an  atmosphere  with  constant  spound  speed,  density  scale  height,  and  Alfvcn 
velocity  (see  McLellan  and  Winterberg  1968  and  NT).  The  plus  modes  all  lie  above  the  upper  dashed  line  12  = A' 
in  Figure  2.  which  corresponds  to  w = c,A  , (the  Lamb  mode'.  The  minus  mode  (there  is  only  one  in  this  case) 
lies  below  the  lower  dashed  line  12  /2A  which  corresponds  to  w = rnkx.  There  are  no  eigenmodes  in  the  region 

between  the  dashed  lines.  This  classification  of  plus  or  minus  modes  refers  here  to  the  character  of  the  eigenmode 


Ll2y2<f>(0)  + (J$Jd:2)\0 
ii22  - A22/1222 


(33) 


F/o.  2.  The  first  five  plus  eigenmodes  and  ttic  only  minus  eigenmode  of  the  two-layer  penumbral  model  for  B-  = 0.5.  y>  = 
y2  --  5/3.  The  quadrangle  represents  the  range  ol  observational  data  (sec  text).  The  crosses  correspond  to  particular  observations 
(Giovanelli  1974). 
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in  the  upper  layer;  in  the  lower  layer,  all  of  the  modes  have  the  character  of  acoustic  waves  (no  magnetic  field  and 
no  buoyancy). 

Observational  data  on  penumbral  waves  are  included  in  Figure  2 for  comparison.  The  most  commonly  reported 
observational  quantities  are  the  period  and  the  horizontal  phase  velocity,  although  they  are  not  always  measured 
simultaneously.  Giovanelli  (1974)  reports  a typical  phase  velocity  of  running  penumbral  waves  of  15  km  s 1 and 
typical  periods  in  the  range  of  180-240  s.  He  did  report  phase  velocities  of  up  to  21  km  s ',  however,  and  gave 
specific  periods  and  wavelengths  for  four  sunspots  (denoted  by  crosses  in  Fig.  2).  The  data  of  Beckers  and  Schultz 
(1972)  appear  to  indicate  a penumbral  oscillation  period  of  255  s.  Moore  and  Tang  (1975)  observed  penumbral 
waves  with  period  270  + 10  s in  a single  sunspot.  Zirin  and  Stein  (1972)  state  that  the  periods  of  penumbral  waves 
in  20  sunspots  were  almost  all  between  240  and  300  s,  and  the  measured  horizontal  phase  velocity  of  9.4  km  s"1 
in  one  spot  was  more  or  less  the  same  in  other  spots  even  when  the  period  varied. 

The  quadrangle  in  Figure  2 represents  the  range  of  observations:  periods  from  180  to  300  s,  and  phase  velocities 
from  9.4  to  21  km  s ',  with  a dashed  line  at  15  km  s 1 to  indicate  the  value  that  Giovanelli  considers  typical. 
The  first  plus  mode  of  the  penumbral  model  passes  through  this  quadrangle.  Although  the  particular  eigenmode 
of  oscillation  of  the  penumbra  is  determined  by  the  excitation,  the  details  of  which  are  uncertain,  the  present 
results  indicate  that  it  is  the  first  plus  eigenmode  that  is  being  excited.  This  agrees  with  our  earlier  conclusions 
(NT). 

The  value /32  = 0.5  used  in  Figure  2 waschosen  to  represent  a typical  penumbra.  In  Figure  3 the  effect  of  changes 
in  £2  on  the  first  plus  mode  is  shown  for  a range  of  fi2  of  two  orders  of  magnitude  (0.05.  0.5.  5.0).  This  constitutes 
a reasonable  set  of  limits  on  )32  for  penumbral  conditions,  and  is  obtained  by  looking  at  the  normal  variation  of 
B (factor  of  4),  ,o  (factor  of  4),  and  c2  (factor  of  1 .5)  expected  in  different  penumbrac.  We  see  that  for  any  reasonable 
value  of  the  penumbral  model  has  a first  plus  mode  within  the  range  of  observations. 

The  vertical  distributions  of  velocity  and  kinetic  energy  of  the  first  plus  mode  (for  jS2  = 0.5)  are  shown  in  Figure 
4 for  A.',  = 0.1995  and  iij  = 0.3628,  corresponding  to  a horizontal  wavelength  A = 3.000  km  and  period  t = 
250  s for  c,2  = 43.5  km2  s'2.  Here  the  nondimensional  height  is  scaled  everywhere  by  the  density  scale  height  in 
the  upper  layer,  //,.  The  velocity  distribution  is  fairly  symmetric,  with  the  maximum  amplitude  occurring  slightly 
above  z = 0 m the  penumbral  photosphere.  The  kinetic  energy,  on  the  other  hand,  is  almost  entirely  trapped  in 
the  lower  layer  (convection  zone)with  maximum  energy  just  below  the  interface.  The  velocity  amplitude  decays 
slowly  with  height  with  a value  of  more  than  25  percent  of  the  maximum  amplitude  at  a distance  of  8 scale  heights 
above  the  level  of  that  maximum. 

There  is  a discrepancy  between  the  height  of  maximum  velocity  predicted  here  (z  T 100  km)  and  that  predicted 
in  our  earlier  paper  (NT.  z ~ 1000  km).  Here  the  Alfven  velocity  increases  exponentially  with  height  above  the 
photosphere,  whereas  in  NT  it  increased  linearly.  Thus,  the  downward  refraction  of  waves  is  much  stronger  in  the 


Fig.  4.  The  distribution  of  vertical  velocity  and  kinetic  energy  of  the  first  plus  cigenniode  with  nondimensional  height  :.H .. 
The  velocity  and  energy  are  each  normalized  to  value  unity  at  the  interface  zjHi  = 0. 
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present  model.  The  actual  situation  is  probably  somewhere  between  these  two  cases.  In  either  case,  the  wave  energy 
lies  mostly  below  the  height  of  maximum  vertical  velocity,  in  the  convection  zone  and  low  photosphere. 


V.  CONCLUSIONS 

The  present  results,  taken  together  with  our  earlier  work  (NT),  indicate  that  running  penumbral  waves  should  be 
identified  with  the  lowest  plus  mode  of  trapped  magneto-atmospheric  waves  in  the  penumbra.  The  vertical  trapping 
is  primarily  due  to  the  increasing  Alfven  velocity  up  into  the  chromosphere  and  the  increasing  sound  speed  down 
into  the  convection  zone.  Most  of  the  energy  of  the  penumbral  waves  lies  in  the  convection  zone  and  low  photo- 
sphere, at  the  same  level  as  the  expected  source  of  excitation  (umbral  oscillatory  convection).  The  maximum  wave 
amplitude  occurs  somewhat  higher. 

The  results  also  indicate  that  penumbral  waves  should  be  observable  in  a photospheric  spectral  line  (see  Fig.  4) 
as  well  as  in  Ha.  There  is  some  indication  of  this  in  the  observation  of  Beckers  and  Schultz  (1972).  Their  data  show 
a 255  s period  oscillation  in  the  penumbra  of  one  sunspot  observed  in  a photospheric  line.  They  present  contours 
of  vertical  velocity  as  a function  of  horizontal  position  and  time  (their  Fig.  1)  in  which  one  may  note  horizontal 
propagation  outward  across  the  penumbra  at  about  the  right  phase  speed.  We  plan  further  observations  in  a search 
for  penumbral  waves  in  the  photosphere. 


Much  of  this  work  was  done  while  we  were  guests  of  the  Max-Planck-Institut  fur  Physik  und  Astrophysik  in 
Mtinchen,  Germany.  We  are  grateful  to  Professor  Biermann  for  the  hospitality  of  the  institut,  and  to  Drs.  H.  U. 
Schmidt,  Friedrich  Meyer,  John  Stewart,  and  Tadashi  Hirayama  for  helpful  discussions.  We  also  thank  Professor 
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ABSTRACT 

Simultaneous  observations  have  been  made  of  velocities  in  the  chromosphere  (in  Hoi  and  in  the 
photosphere  (in  the  nonmagnetic  Ee  l line  X5576)  of  three  sunspots.  The  results  reveal  waves  propa- 
gating horizontally  outward  across  the  penumbra  in  the  photosphere  with  about  the  same  period 
as  the  running  penumbral  waves  in  Ho  (250  290  s).  The  photospheric  waives  are  more  intermittent 
and  have  higher  horizontal  phase  velocity  (by  a factor  of  2 or  more)  than  the  chromospheric  penum- 
bral waves. 

Subject  headings:  hydromagnetics  — Sun:  atmospheric  motions  — - Sun:  sunspots 


t.  INTRODUCTION 

One  of  the  most  interesting  recent  discoveries  con 
corning  the  dynamics  of  the  solar  atmosphere  is  that  of 
waves  observed  in  Ha  that  propagate  outward  across 
sunspot  penumbrae.  These  “running  penumbral  waves” 
(or  simply  “penumbral  waves”)  were  discovered  inde- 
pendently by  Zirin  and  Stein  (1972)  and  by  Giovanelli 
( 1972,  1974).  They  tire  often  seen  in  large  sunspots  with 
a stable,  regular  penumbra,  and  their  observed  hori- 
zontal phase  velocities  are  in  the  range  8 -25  km  s"1, 
periods  in  the  range  150-270  s,  and  horizontal  wave- 
lengths in  the  range  2500-5800  km.  The  associated 
motions  are  largely  vertical  ( i.e. , normal  to  the  solar 
surface). 

The  nature  of  the  wave  motion  indicates  that  com- 
pressibility, buoyancy,  and  magnetic  field  all  play  a role 
in  penumbral  waves.  The  observed  phase  velocities  are 
indicative  of  conditions  in  the  penumbral  photosphere 
rather  than  the  chromosphere,  t Eor  example,  the  Alfven 
speed  at  height  of  formation  of  Ho  is  much  larger  than 
25  km  s E.  This  suggests  that  penumbral  waves  are  a 
photospheric  phenomenon.  We  have  given  a theoretical 
description  (N ye  and  Thomas  1974,  1<>70  > of  penumbral 
waves  as  magnetoatmosphcric  waves  that  are  vertically 
trapped  at  photospheric  levels,  the  trapping  being  due 
to  refrac  tion  by  the  increasing  Alfven  speed  up  into  the 
chromosphere  and  the  increasing  sound  speed  down 
into  the  convection  zone.  The  waves  are  evanescent  at 
heights  of  forma. ion  of  Ho,  but  have  their  largest  ampli- 
tude there  due  to  the  rapid  decrease  in  density  with 
height.  The  wave  energy  is  mostly  in  the  photosphere 
and  subphotosphere,  which  is  also  the  location  of  the 
probable  source  of  excitation  of  the  waves,  nnmclv, 
overstable  conviction  in  the  umbra  (Musman  1967; 
Sit vage  1969;  Moore  1975) 
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For  a maximum  vertical  velocity  of  order  1 km  s 1 
in  Ha  (Giovanelli  1974),  the  theory  (Nye  anil  Thomas 
1974,  1976)  predicts  vertical  velocities  of  order  100  m s 1 
or  greater  in  the  photosphere;  thus,  penumbral  waves 
should  also  be  observable  in  a pho'ospherii  spectral  line. 
This  prediction  prompted  the  simultaneous  photo 
spheric  and  chromospheric  observations  reported  here. 

II.  OBSERVATIONS 

We  used  the  Sacramento  1’eak  vacuum  tower  tele- 
scope (Dunn  1969),  echelle  spectrograph  (Dunn  1971 
and  the  diode  array  (Dunn  and  Spence  1975  to  obtain 
spectral  scans  in  a manner  similar  to  that  of  Musman 
(1974).  At  each  position  of  a two-dimensional  raster  w< 
recorded  on  magnetic  tape  the  intensity  measured  at 
4(1  equally  spaced  wavelengths  in  t spectral  interval  of 
1 A centered  about  the  nonmagnetic  line  of  Ee  t at 
X5576.099.  This  line  is  formed  in  the  upper  penumbral 
photosphere.  Alt  rock  cl  al.  (1975)  s ate  that  in  the  quiet 
photosphere  90  percent  of  the  lint  center  contribution 
function  is  formed  between  250  and  550  km  above 
optical  depth  unity  at  5000  A.  We  assume  a similar 
behavior  for  the  penumbra.  Each  measurement  is  lint 
ited  by  a 1"  (700  km)  aperture  and  s performed  in  2 15 
of  a second.  Each  scan  line  contains  40  positions  with 
1"  spacing.  Each  raster  contains  ,'nrcc  adjacent  lines 
with  1"  spacing  and  takes  25.4  r to  complete.  The 
raster  was  positioned  visually  such  that  the  center  scan 
line  passed  through  the  center  of  the  sunspot  umbra  in 
each  observation;  in  this  way,  propagation  along  the 
scan  lines  corresponds  to  radial  propagation  across  the 
sunspot.  A ty  pical  run  consisted  of  64  scans  covering  a 
time  period  of  27  minutes.  After  each  scan  the  image  is 
returned  to  the  center  of  the  raster  and  an  Ho  photo- 
graph of  the  slit  jaws  is  'aken.  Th  s process  occupies 
about  5 s of  the  total  scan  time.  The w Ho  photographs 
were  later  combined  to  form  a movie. 

Velocities  are  obtained  by  locating  the  center  of  the 
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line  profile  by  interpolation  on  a smoothed  line  profile 
as  described  by  Musnian  (1974).  All  velocities  are  re 
ferred  to  the  mean  velocity  of  each  raster.  We  tveraged 
the  velocities  of  the  three  adjacent  lines  in  order  to  sup- 
press the  noise  in  a single  measurement  and  obtained 
the  equivalent  of  one  40-point  scan  line.  Since  the  veioe 
it \ is  measured  at  a slightly  different  time  and  we  scan 
an  odd  number  of  lines  in  alternate  directions,  a per- 
fectly uniform  oscillatory  pattern  on  the  Sun  c mid  give 
rise  to  an  apparent  phase  velocity  in  the  direction  of 
scanning  of  the  odd-numbered  lines.  However  this  ap- 
parent phase  velocity  is  greater  than  3(XX)  krt  s 1 and 
introduces  a negligible  error  in  measuring  phase  veloc- 
ities two  orders  of  magnitude  smaller. 

We  observed  two  sunspots  in  McMath  regions  13875 
and  13890  on  1975  October  5 9 and  October  16  19.  I’re 
liminary  observations  were  also  obtained  for  a sunspot 
in  McMath  region  13738  on  1975  June  28- July  1,  but 
detailed  results  will  be  given  here  only  for  the  later  two 
spots. 

III.  RESULTS 
<j)  Chromosphere 

When  the  Ha  slit-jaw  pictures  are  projected  as  a 
movie,  the  running  penumbral  waves  are  most  con- 
spicuous for  McMath  region  13875.  The  waves  appear 
continuously  for  all  observations  of  this  sunspot  over  a 
period  of  5 days.  They  have  an  average  period  of  240  ± 
10  s and  propagation  velocities  in  the  range  from  10- 
18  km  s*1.  These  waves  then  appear  to  be  typical,  con- 
sistent with  the  observations  of  Zirin  and  Stein  (1972) 
and  Giovanelli  (1972,  1974). 

On  viewing  the  Ha  movie  of  McMath  region  13890, 
one  has  the  impression  of  waves  propagating  outward 
across  the  penumbra,  but  the  pattern  is  not  distinct 
enough  to  make  measurements  of  periods  or  velocities. 
Thtre  is  no  indication  of  running  penumbral  waves  in 
McMath  region  13738,  but  the  time  interval  between 
pictures  is  longer  in  that  observation,  making  the  waves 
more  difficult  to  detect. 

ft)  Photosphere 

To  determine  the  reliability  of  the  X5576  velocity 
data,  it  was  first  necessary  to  determine  the  noise  in  the 
signal.  To  do  this,  a series  of  scans  was  made  out  of 
focus,  giving  a uniform  solar  signal.  From  this  the  noise 
was  found  to  have  nearly  equal  power  at  all  frequencies 
with  rms  velocity  on  the  order  of  26  ms*1.  We  also 
measured  a nearly  uniform  spectrograph  drift  of  21  mA 
hr  *.  The  average  value  of  each  scan  was  subtracted  to 
remove  the  effect  of  the  spectrograph  drift.  Next,  the 
time  series  for  each  of  the  40  points  in  the  scan  was 
analyzed.  The  mean  for  each  time  series  was  removed. 
This  eliminated  the  i mspicuous  effect  of  Evershed  flow. 
The  best  estimate  of  the  true  signal  was  obtained  by 
using  a bandpass  fillet  in  the  manner  described  by 
Hrault  and  White  (1971  . The  filter  removed  high-fre- 
quency (period  < 70s)  noise  plus  a very  low-frequency 
(period  > 800  s)  drib  due  to  an  instrumental  problem 
in  the  diode  array  The  power  spectra  for  the  unfiltered 
data  are  essentially  Hat  for  high  frequencies,  which  is 


indicative  of  white  noise  and  consistent  with  the  out- 
of-focus  data. 

The  filtered  data  clearly  show  5 minute  oscillations  in 
tne  portion  of  the  scan  outside  of  the  spot.  The  data 
also  show  5-minute  oscillations  in  the  umbra,  with 
typical  amplitudes  of  50  ms  but  since  these  are 
usually  180°  out  of  phase  with  the  oscillations  outside 
the  spot,  we  suspect  that  this  is  a systematic  error  due 
to  referring  all  velocities  to  the  average  of  each  scan. 
Repeating  one  of  the  scans  with  a velocity  referenced 
to  the  average  umbra!  velocity  in  an  attempt  to  remove 
this  effect  does  not  significantly  change  the  results  re- 
ported below.  The  error  in  velocity  reference  is  opposite 
to  the  effect  of  scattered  light,  which  would  give  oscilla- 
tions in  phase.  We  looked  primarily  for  horizontally 
propagating  waves  in  the  penumbra  resembling  those 
observed  in  Ha.  Neither  scattered  light  nor  an  inaccu- 
rate velocity  reference  could  introduce  such  a systematic 
propagation. 

The  filtered  data  show  oscillations  in  the  penumbral 
photosphere,  but  not  in  the  form  of  a continuous  out- 
ward-propagating wave  train  sometimes  seen  in  Ha. 
There  do  appear,  on  four  occasions,  what  we  term 
penumbral  photospheric  events.  These  are  isolated  wave 
packets  that  propagate  outward  across  the  penumbra. 
The  packets  consist  of  one  or  two  complete  oscillations 
whose  period  (250  290  s)  is  comparable  to  that  of  the 
chromospheric  penumbral  waves.  Figure  1 shows  an 
example  of  one  such  event,  as  illustrated  by  the  time 
history  of  the  line-of-sight  velocity  at  each  scan  point 
across  the  penumbra.  The  heliocentric  angle  of  the  sun- 


Kig.  1 .— - An  example  of  a penumbral  photospheric  event  (event 
3 in  Table  1).  The  time  records  of  the  filtered  velocity  signal  at 
points  across  the  penumbra  are  shown.  The  numbered  scan  posi- 
tions (with  1"  spacing)  correspond  to  those  shown  in  Fig.  2. 
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spot  during  each  of  the  four  events  was  roughly  .<0°,  so 
that  line-of-sight  velocities  were  predominantly  vertical. 
The  amplitude  of  each  of  the  four  events  is  nearly  an 
order  of  magnitude  larger  than  the  rms  noise.  There 
may  have  been  a fifth  event  in  McMath  region  1.1738, 
but  it  occurred  too  near  the  beginning  of  the  record  to 
be  identified  unambiguously. 

The  horizontal  phase  velocities  of  the  four  events  are 
in  the  range  17  -90  km  s *.  Individual  phase  velocities, 
amplitudes,  periods,  locations,  dates,  and  times  of  the 
four  events  are  listed  in  Table  1.  The  locations  given  in 
Table  1 refer  to  the  numbered  positions  across  the  sun 
spot  shown  in  Figure  2,  which  also  shows  plots  of  the 
continuum  intensity  across  the  scan.  The  times  given 
in  Table  1 refer  to  the  central  time  of  each  event.  This 
is  the  time  when  the  maximum  amplitude  of  each  wave- 
packet  passed  the  central  location  in  its  range  of  propa- 
gation. The  periods  listed  in  Table  1 are  consistent  with 
the  data  of  Beckers  and  Schultz  (1172),  who  found 
power  at  255  s in  the  penumbral  photosphere  of  the  sun- 
spot they  studied. 

To  determine  whether  these  penumbra!  photospheric 
events  could  have  been  due  to  chance,  we  calculated 


McMath  Region  13890 


the  probability  that  the  phase  of  a sequence  of  inde- 
pendent oscillators  with  the  same  frequency  but  ran- 
dom phase  could  combine  to  form  an  apparently  propa- 
gating wave.  This  type  of  problem  has  been  treated  bv 
Roach  (1968).  The  probability  of  having  the  phase 
differences  of  five  random  oscillators  out  of  40  within 
20”  is  5 X 10  4.  We  analyzed  in  detail  six  independent 
sets  of  data  from  six  different  days  and  two  sunspots, 
for  a total  time  of  2.7  hours.  Al!  four  events  extended 
over  live  to  seven  positions  and  are  statistically  sig- 
nificant. 

IV.  CONCLUSIONS 

Our  observations  at  chromospheric  and  photospheric 
levels  in  three  sunspots  lead  us  to  the  following  conclu- 
sions: 

1.  Typical  running  penumbral  waves  were  observed 
at  chromospheric  levels. 

2.  At  photospheric  levels  we  observed  statistically 
significant  events  consisting  of  isolated  wave  packets 
(one  or  two  complete  oscillations)  propagating  outward 
across  the  penumbra.  Four  events  were  seen  during  2.7 
hours  of  observation.  The  period  of  these  oscillations  is 
comparable  to  that  of  the  penumbral  waves  in  Ha.  The 
existence  of  these  events  is  a partial  confirmation  of  the 
prediction  of  Nye  and  Thomas  (1974,  1976). 

3.  The  waves  observed  in  the  penumbral  photosphere 
are  more  intermittent  and  have  higher  horizontal  phase 
velocity  (by  a factor  of  2 or  3 or  more)  than  the  penum 
bral  waves  in  Ha.  Thus,  the  connection  (if  any)  between 
the  penumbral  photospheric  events  and  the  chromo- 
spheric penumbral  waves  is  unclear  at  present,  and  in 
any  event  is  more  complicated  than  the  resonant  mode 
model  of  Nye  and  Thomas  (1976). 


McMath  Region  13875 

Fig.  2.— Continuum  intensity  versus  position  across  the  scan 
for  the  two  sunspots  in  Table  1.  The  numbered  scan  positions  (0- 
40)  have  1"  spacing. 
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TABLE  1 

Parameters  for  the  Photospheric  Penumbral  Events 


Horizontal 

Phase 


Sunspot 

Location* 

Velocity 
(km  s“‘) 

Amplitude 

(m  s-‘) 

Period 

(s) 

Date  and 
Time  ( CT  • + 

13875 

21  -27 

<81 

220 

204 

10/6  75, 
14:27:52 

13875 

11-7 

04 

255 

258 

10/6/75, 

14:22:08 

13890 

25-31 

37.5 

265 

252 

10/16/75, 
14:42 :54 

13890 

18-13 

47 

170 

290 

10/16/75, 

14:26:53 

* Location  refers  to  the  numbered  scan  positions  in  Fig.  2. 
t Time  refers  to  the  central  time  of  the  event. 


I 


40 


1.178  M USMAN,  NVE,  AND  THOMAS 


KEIT'.KEXCES 

Alt  rock,  R.  November,  L.  J . Simon,  (J.  W.,  Nlilkey,  R.  \\ 

and  Worden,  S.  IV  1075,  Solar  Tliys.,  43,  33. 

Beckers,  J.  M.,  and  Schultz,  R.  B.  1072,  Solar  Tliys.,  27,  61. 

Brault,  J.  W.,  and  White,  ().  R.  1071,  A sir.  and  Ap.,  13,  160 
Dunn,  R.  B 1060,  Sky  and  Tel.,  38,  368. 

1971,  in  The  Stenzel  Symposium  on  Solar  Physics,  Atomic 
Spectra,  and  ( iaseous  Xebulae,  ed.  K.  B.  Gehhic  (NBS  Spec. 

Bul»  No.  353),  p.  71 . 

Dunn,  K.  B . and  Spence,  <i.  K.  1073,  Bull  A AS,  5,  271. 

Ciiovanelli,  R.  G.  1072,  Solar  Tliys.,  27,  71. 

. 1074,  in  I AT  Symposium  Mo.  50,  Chromospheric  Tine 

Structure,  ed.  R G.  Athay  (Dordrecht:  Keidel),  p.  137. 

Steven  Musman  and  Alan  H.  Nve:  Sacramento  Peak  Observatory,  Sunspot,  NM  88349 

John  \.  Thomas:  Department  of  Mechanical  and  Aerospace  Sciences,  University  of  Rochester,  Rochester,  NY  14627 


Moore,  R.  L.  1073,  Solar  Tliys.,  30,  403. 

Musman,  S.  1067,  Ap.  J .,  149,  201. 

. 1074,  Solar  Tliys.,  36,  313. 

Nve,  A.  H.,  and  Thomas,  J.  H.  1074,  Solar  Tliys.,  38,  300. 

. 1076,  Ap.  J.,  204,  582. 

Roach,  S.  A.  1068,  The  Theory  of  Random  Clumping  (London: 
Methuen),  p.  6. 

Savage,  B D.  1060,  A p.  J .,  156,  707. 

Zirin,  1L,  and  Stein,  A.  1072,  Ap.  J.  (letters),  178,  L85. 


I 


SOLAR  MAGNETO-ATMOSPHERIC  WAVES 

John  H.  Thomas 

Department  of  Mechanical  and  Aerospace  Sciences 

I and 

C.E.  Kenneth  Mees  Observatory 

i 

University  of  Rochester 
Rochester,  New  York  14627,  U.S.A. 

(Paper  to  be  presented  at  IAU  Colloquium  36,  "Energy  Balance  and  Hydro- 
dynamics of  the  Chromosphere  and  Corona,"  Nice,  France,  September  6-10,  1976) 

***** 

I would  like  to  say  a few  words  about  the  theory  of  what  I call 
magneto-atmospheric  waves,  which  also  go  under  the  name  of  magneto- 
acoustic-gravity waves.  These  are  waves  that  involve  the  combined  restoring 
forces  due  to  compressibility,  gravity,  and  the  magnetic  field.  For  many 
waves  observed  in  the  solar  atmosphere,  all  three  of  these  effects  are 

important 

Nearly  all  work  on  magneto-atmospheric  waves  has  been  based  on 
linearized  theory  for  a plane-stratified,  nondissipative  atmosphere. 

(See  Chiu  1971  for  a discussion  of  the  nonlinear  case.)  Wave  propagation 
is  then  determined  by  the  vertical  distribution  of  three  parameters: 
the  sound  speed  c,  the  AlfvSn  speed  vA,  and  the  local  density  scale 
height  H.  The  theory  is  complicated  by  the  fact  that  the  medium  is  both 
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anisotropic  and  inhomogeneous.  Gravity  and  the  magnetic  field  each 
Introduce  a preferred  direction.  The  inhomogeneity  Influences  wave 
propagation  most  strongly  through  the  rapid  increase  of  the  AlfvGn 
speed  with  height. 

The  simplest  approach  to  atmospheric  wave  problems  is  to  find 
a dispersion  relation  for  plane  waves.  However,  this  approach  has 
not  been  very  fruitful  for  solar  magneto-atmospheric  waves.  Unless 
all  three  wave  parameters  c,  u^,  and  H are  constant  with  height,  we 
can  find  only  a "local"  dispersion  relation  (MacDonald  1961;  McLellan 
and  Winterberg  1968;  see  also  Bel  and  Mein  1971;  Michalitsanos  1973; 
Nakagawa  et  al_.  1973;  Yeh  1974)  that  holds  only  for  waves  whose  vertical 
wavelsng-.h  is  much  smaller  than  the  smallest  scale  height  of  the  un- 
disturbed atmosphere.  We  can  see  how  restrictive  this  approach  is  by 
recalling  that  in  the  photosphere  and  much  of  the  chromosphere  the 
density  scale  height  is  smaller  than  our  limit  of  observational 
resolution. 

There  is  one  special  case,  first  studied  by  C.P.  Yu  (1965;  see  also 
Chen  and  Lykoudis  1972;  Nye  and  Thomas  1974),  in  which  a,  v^,  and  H are 
all  constant.  This  is  the  case  of  an  isothermal  atmosphere  with  a 
horizontal  magnetic  field  that  decreases  exponentially  with  height  as 
exp  ( -z/2H ).  This  leads  to  a global  dispersion  relation.  However, 
magnetic  fields  on  the  sun  usually  vary  more  slowly  than  the  density, 
so  this  case  has  limited  applicability. 

Perhaps  the  simplest  model  atmosphere  for  solar  magneto-atmospheric 
waves  is  an  isothermal  atmosphere  with  a uniform  magnetic  field.  In 
this  case  the  sound  speed  a and  the  density  scale  height  H are  constant. 
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This 


but  the  Alfv§n  speed  increases  exponentially  with  height  as  eZ//2}1 . 
rapid  increase  of  AlfvGn  speed  with  height  causes  strong  refraction  and 
downward  reflection  of  waves,  which  can  not  be  adequately  discussed  with 
the  local  dispersion  relation.  We  must  go  back  to  the  basic  wave 
equations  with  variable  coefficients,  for  which  solutions  have  been 
found  in  certain  cases. 

Approximate  solutions  for  the  case  of  a uniform  vertical  magnetic 
field  were  first  given  by  Ferraro  and  Plumpton  (1958;  see  also  Weymann 
and  Howard  1958;  Stepien  1967).  There  are  three  types  of  wave  modes 
in  this  case.  One  is  an  AlfvSn  wave  having  only  nondivergent  horizontal 
motions.  The  other  two  modes  have  both  vertical  and  horizontal  motions, 
and  may  be  described  as  gravity-modified  magneto-acoustic  waves,  or  as 
magnetically-modified  acoustic-gravity  waves,  depending  on  the  relative 
valves  of  the  sound  speed  and  AlfvSn  speed  at  some  reference  level. 

Meyer  (1968)  used  an  approximate  solution  similar  to  Ferraro  and  Plumpton's 
in  a model  for  flare-induced  coronal  waves,  or  Moreton  waves.  In  his  model 
the  chromosphere-corona  transition  region  provides  a lower  reflective 
boundary,  and  the  resulting  trapped  modes  of  gravity-modified  fast  waves 
are  identified  with  the  coronal  waves. 

For  the  case  of  an  isothermal  atmosphere  with  a uniform  horizontal 
magnetic  field,  Alan  Nye  and  I (Nye  and  Thomas  1976a)  have  recently 
given  an  e;;act  solution  for  the  wave  modes.  The  wave  equation  in  this 
case  transforms  to  the  hypergeometric  equation,  and  solutions  can  be 
expressed  in  terms  of  hypergeometric  functions.  We  applied  this 
solution  to  a model  of  flare-induced  coronal  waves  similar  to  Meyer's, 


but  with  the  magnetic  field  horizontal  instead  of  vertical.  The 
resulting  trapped  modes  propagate  horizontally,  with  very  little 
dispersion  for  all  horizontal  wavelengths  less  than  about  100,000  km 
(see  Fig.  1).  Thus  a flare-induced  pulse  will  propagate  horizontally 
over  large  distances  with  little  distortion,  as  is  indeed  observed. 

We  have  also  used  our  exact  solution  in  a model  for  running 
penumbral  waves  (Nye  and  Thomas  1976b).  The  model  consists  of  an 
upper  isothermal  layer  with  a uniform  horizontal  magnetic  field, 
representing  the  penumbral  photosphere  and  low  chromosphere,  and  a 
lower  aciabatic  layer  with  no  magnetic  field,  representing  the  under- 
lying convection  zone.  The  penumbral  waves  are  identified  with  fast 
modes  that  are  trapped  by  the  increasing  AlfvGn  speed  above  and  the 
increasing  sound  speed  in  the  convection  zone  below.  Figure  2 is 
a diagnostic  diagram  showing  computed  eigenmodes  for  a typical 

penumbra,  along  with  the  range  of  observed  frequency  and  wavelength 
of  penumbral  waves.  There  is  good  agreement  between  observations  and 
the  lowest  eigenmode.  Based  on  the  form  of  the  eigenmode,  the  theory 
predicts  that  penumbral  waves,  which  were  discovered  in  Ha,  should  also 
be  observable  in  the  photosphere.  This  prediction  was  recently  confirmed 
by  Steve  Musman,  Alan  Nye,  and  myself  (Musman,  Nye,  and  Thomas  1976)  in 
observations  with  the  tower  telescope  and  diode  array  at  Sacramento 
Peak  Observatory. 

I want  to  mention  briefly  one  other  aspect  of  magneto-atmospheric 
waves  which  I have  been  working  on  recently.  In  the  case  of  a horizontal 
magnetic  field  there  is  a singularity  in  the  basic  wave  equation  at  a 


height  at  which 
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When  a or  vA,  or  both,  vary  with  height,  then  for  slow  waves  of  a 
certain  frequency  w and  horizontal  wavenumber  k there  exists  a 
critical  height  zc  at  which  this  relation  is  satisfied.  In  the  local, 
or  WKB  approximation,  the  vertical  wavenumber  becomes  infinite  at 
the  critical  level,  and  waves  take  an  infinite  time  to  reach  this  level; 
thus,  the  critical  level  acts  as  an  absorbing  barrierfor  waves.  But  in 
fact  the  WKB  approximation  breaks  down  at  the  critical  level,  so  this  is 
not  an  adequate  discussion. 

This  phenomenon  is  similar  to  the  critical  levels  that  occur  in 
the  stability  of  inviscid  shear  flow  (Lin  1966)  and  in  the  propagation 
of  internal  gravity  waves  in  a shear  flow  (Booker  and  Bretherton  1967). 

In  the  present  case  the  singularity  is  logarithmic  and  the  singular 
point  is  thus  a branch  point.  We  can  use  the  technique  of  Booker  and 
Bretherton  to  decide  on  the  proper  branch  of  the  solutions  and  to  connect 
solutions  above  and  below  the  critical  level.  The  conclusion  is  that 


waves  are  not  completely  absorbed,  but  are  attenuated  as  they  pass  through 
the  critical  level.  The  attenuation  is  algebraic,  rather  than  exponential 
as  in  the  case  of  internal  gravity  waves  in  a shear  flow.  Energy  is 
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Fig.  I— Phase  velocity  v,  and  group  velocity  u,  of  the  first  mode  of  coronal  oscillation  plotted  as  a function  of  horizontal 
wavenumber. 
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Fir.,  2.  The  first  five  plus  eigcnmodes  and  the  only  minus  cigcnmodc  of  the  two-layer  pcnumbral  model  for  P ‘ = 0 5,  y,  = 
y:  = 5/J.  The  quadrangle  represents  she  range  of  observational  data  (sec  test)  The  crosses  correspond  to  particular  observations 
(Uiovanclli  1974). 
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